Silllif 










: Wi 




NMMMil 

Qmmmmm i 

V*T-- ; S 


. ^ 


I! 


I! 


i^AC^^QN NUMBER) O 'U § -*■ 

(THRU) 

i£fi 

1 


(PAGES) 


IMZIM 

(NASA CR OR TMX OR AD NUJ 


NUMBER) 


(CODE) 

2^3 


(CATEGORY) 




■ 




mm - fill ilillili 

; .■ . : 









28 March 1968 


RSIC-769 


RECENT ADVANCES IN COHERENT OPTICS: 
FILTERING OF SPATIAL FREQUENCIES; HOLOGRAPHY 

by 

Serge Lowenthal 
Yves Belvaux 


Revue d’Optique, Theorique et Instrumentale, 

46, No. 1, pp. 1-64(1967) 


Translated from the French 


THIS DOCUMENT HAS BEEN APPROVED FOR PUBLIC RELEASE 
AND SALE; ITS DISTRIBUTION IS UNLIMITED. 


Translation Branch 

Redstone Scientific Information Center 
Research and Development Directorate 
U. S. Army Missile Command 
Redstone Arsenal, Alabama 35809 



RECENT ADVANCES IN COHERENT OPTICS: 
FILTERING OF SPATIAL FREQUENCIES; HOLOGRAPHY 

by 

Serge Lowenthal 
Yves Belvaux 


Compagnie Generate de Telegraphic Sans Fil (CSF) 
(General Wireless Telegraphy Company) 


New applications of coherent light are presented: 

a) Holography: recording of three-dimensional objects; diffuse illumi- 
nation interferometry; reproduction of holograms; Fourier-transform holography 
and geometrical properties of holograms; 

b) Spatial filtering: optical data processing, especially pattern recogni- 
tion; automatic reading of characters, identification of fingerprints and other 
signals. 


INTRODUCTION 

Coherent optics is not a new field of physics: ever since Fresnel at the 
beginning of the 19th century there has been a succession of studies on luminous 
interference and diffraction phenomena. 

However, the appearance in 1961 of the gas laser marks a discontinuity 
in this movement because of the rapid and considerable progress which has been 
accomplished since then in the field of coherent optics. 

Two methods have particularly benefited from the increase in luminance, 
and in temporal and spatial coherence brought about by laser, compared with 
the classical sources: 

a) The filtering of spatial frequencies in coherent light, thanks to the 
work of A. Marechal (1953) [1] ; 

b) Holography, conceived by D. Gabor in 1948 [2 ] and subsequently 
perfected by the researchers of the University of Michigan [3,4], 
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In this paper we shall discuss the results of recent studies carried out 
under a DEME contract (Note: Contract 65-34-276 entitled "The Possibilities 
of Filtering of Spatial Frequencies") whose aim was to explore certain current 
possibilities of these methods. 

In Chapter I we shall recapitulate certain elementary concepts relating 
to diffraction, in a form adapted to the problems treated subsequently. 

Chapter II is devoted to the holograms. We know that the latter permit 
recording and then restitution on a photographic emulsion, with the aid of a 
coherent background, both the phase and the amplitude of a coherent light wave. 

The results which we have obtained relate to the three-dimensional 
holograms (images in relief) , to the duplication of holograms, to interferometry 
in diffuse light, to the formation of images by holography, and to Fourier- 
transform holography. 

In Chapters III and IV we shall first recall die principles of the filtering 
of spatial frequencies, and will show that this method, based on diffraction in 
the optical instruments — a method which has yielded spectacular results in the 
improvement of optical images [ 5-7 ] — is in fact a method of data processing. 

In effect, the coherent optical filtering implies a convolution relationship 
between two functions . 

Convolution is a linear functional which comprises some important, 
special cases, notably the correlation, integration, differentiation, and scalar 
product of two functions. 

It can be seen that by means of this mathematical apparatus it is possible 
to carry out data processing. It is this aspect of the problem which we have 
specially developed herein. 

The fact that, on the one hand, this treatment is carried out in two 
dimensions (x,y) and that, on the other hand, it requires only a static filtering 
by a material screen whose transparency is suitably modulated, endows this 
method with great simplicity. Thus the optical correlation of two radar charts 
is as easy as that of two rectangular gating pulses. 

In this way coherent optics makes its entry into an area which up to now 
has been reserved for electronic computers. 

In Chapter IV the optical filtering is applied to pattern recognition. 

This term actually comprises some highly diverse processes, all of which, 
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however, are treated by an analogous formalism. A satellite detection system, 
for example, gives a signal from which a number of data are expected such as 
presence, position or velocity; this is a form of pattern recognition. 

Among the results presented herein, we have stressed the automatic 
recognition of characters in a printed text. In effect, despite its apparent 
simplicity, this problem is beset with a number of traps whose solution makes 
it possible to approach other problems, such as the treatment of electric signals 
after spatial transposition, with a greater hope of success. 

The study presented herein makes no claims for being exhaustive; the 
subject is in the process of development, as is shown by the bibliographical 
references given at the end of this article. 

I. GENERAL DISCUSSION OF DIFFRACTION 


1. Introduction 

Coherent optics is based on the phenomena of diffraction. Below we 
shall discuss some of the elements of the theory of diffraction used in this study, 
particularly the form of the correspondence between two planes, one of them 
comprising an illuminated aperture which diffracts the light, and the other being 
the observation plane. 

In the mathematical expressions for this correspondence there appear 
certain pure phase factors which are frequently neglected, because they dis- 
appear during recording done by means of detectors such as the photographic 
plate. 

In holography and in the filtering of spatial frequencies these terms 
play an important role: it is for this reason that an effort is made to show 
them in the formulas which will be developed in the case of 

a) Fraunhofer's diffraction, relating to observation in the vicinity of 
a focus, or at infinity; 

b) Fresnel's diffraction, relating to all other cases. 

These two cases are deduced from the analytical expression of Huygens' 
principle. 
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2. The Huygens-Fresne! Transformation Principle 

This principle may be stated as follows: 


If the electromagnetic field (produced by a coherent source) is known 
on a closed surface S which does not contain any sources in its interior 
(Figure 1) , it is known at any point P in the interior of surface S. 



Figure 1. Huygens’ Principle. If the 
electromagnetic field is known on a 
closed surface S not containing any 
sources such as 2 , it is known at any 
point P located within S. 


It is in this form that Huygens' 
principle will be used in the overall study 
of the holograms. 

The analytical expression of the 
phenomena is easy only when the surface 
is an absorbing plane n o pierced by an 

aperture. Plane n o is completed by the 

half-sphere of infinity whose contribution 
at P is zero. 

For an optical instrument the dif- 
fracting opening is the exit diaphragm. 


Let n be the diffracting plane 

( Figure 2) . We shall now determine the 
complex amplitude of the electromagnetic field diffracted by n at point P 
situated in a plane parallel to n . 



Figure 2. Fresnel’s Diffraction. Il 0 
is the diffracting plane; is the obser- 
vation plane. 


This result can be obtained by means of several expressions. We shall 
use Sommerfeld's formula [8 ] : 
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(1.1) g(x,y,z) = 


e C C uf \ ikr Ak l\, , 

— I I hfx ,y )e cos 01 — - -5 jdx dy , 

2ir J J n Q \ .0 0 / \r r J o J o 

where h , y^ j is the complex amplitude at a point M q having coordinates X q , y Q 

in plane II q of origin O; g(x,y, z) is the diffracted amplitude at point P(x,y,z) , 

r is the distance M P, 0 is the angle between M P and the normal, M N to n , 
o 0 o o 

and k = 2ir/\, the wave number . 

Formula (1. 1) , valid in the scalar approximation for a quasi- 
monochromatic wave of pulsation oj , is simplified when the diffraction angles are 
small. We then have cos 9 - 1 . If, moreover, the distance z = d between 
and n Q is large compared with X, we can neglect 1/r 2 compared with ik/r , and 

take this last term outside the integral, because r — z = d, whence, we have 

( 1 . 2 ) g(x,y.d) = h(x o .y o ')e lkr dx o dy o . 

o 


ikp 

In factor e , we have 


(1.3) 


r = 

d 2 + ( x - x 

2 + ( 

y - y V 


V 0 / 

\ 

0/ 


V. 


Now, d is large compared with the other terms, and (1.3) may be expanded: 

1 


(1.4) 


d + 


2 d 


( x - X o) ! + ( y - y o) ! 


and expression (1.2) becomes 
(1.5) 


g ( x, y , d) = -^o lkd // IIo h (v y o) 



fik 


2 , / \ 2 1 

x exp' 

I2d 

(* - X o) 

+ ( y “ y o)J 


►dx dy . 
o o 


This expression will serve as the point of departure for establishing the 
properties of both Fresnel's and Fraunhofer's diffraction. 
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Vectorial Notation. Vector OM of components x , y , contained in 
o o o 

plane n , will be denoted as M . Likewise vector SP of components x,y, con- 
o o 

tained in plane T will be denoted as P. 

We shall also write 

h ( X o ,y o) = h ( M o)’ g(x,y,d) = g(P,d) , 

dx dy = dM , dxdy = dP. 

oo o 

With these conventions, the term ^x - x^ 2 + ^y - y^ 2 becomes the scalar 

square ^P - M^ 2 and when expanded becomes P 2 + M^ 2 - 2PM q , the term 

PM being the scalar product of P and M . 
o o 

The analytical form (I. 5) of the Huygens- Fresnel principle is therefore 
written as 

dM . 
o 


(1.6) g(P,d> = -ie ikd / n »(M o )exp 2j (P - M g ) 


3. Fresnel Diffraction and Fresnel Transformation 

Expression (1. 6) gives the complex amplitude of the field diffracted 
by plane II at a point P situated at a finite distance from n^; hence this is a 

Fresnel phenomenon. 


When the observation is made in a plane parallel to n^, the factor in 

front of the integral is a constant which may be neglected; hence, by writing 
g(P) forg(P,d), (1.6) becomes 


(1.7) 


s( p > = /„ h(M o )e X p | (P - M o J 


dM 


which can be further written as 
ikP 2 / 2d 


(1.8) g(P) = e 


4 h ( M o) exp i( M o ! - 2P - M o) 


dM . 
o 
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If phenomenon (I. 8) is recorded by photographic plate which is sensitive only 
to illumination lg(P) I 2 , the quadratic factor exp (ikP z /2d) disappears and may 
be neglected in the formulas. By contrast, in holography this factor is recorded, 
and hence the complete expression (1. 7) must be preserved. 

Relation (1.7) is a convolution product symbolically denoted as 

(1.9) g(P) = h(P) * e lkI>2/2d . 

This expression is a Fresnel transformation. It can be shown that the 
latter is reciprocal, i. e. , if we have (I. 9) we also have 

(L10) h ( M o) = W g ( M o) ¥e ° 

The two relations (1.9) and (1. 10) express, for the Fresnel diffraction, 
the correspondence between a diffracting plane n and an observation plane JF; 
the roles of the two planes may be reversed. 


Special Case. The diffracting plane is reduced to a luminous point M, 
and then h^M^ is represented by the Dirac distribution 6^M - and the 
distribution of the complex luminous amplitudes on (F becomes 


( 1 . 11 ) 


g(P) = exp 


£< m - p > ! 
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a result which is readily obtained also by a direct calculation. 


4. Fraunhofer's Diffraction and Fourier Transformation 

Fraunhofer's diffraction takes place in the vicinity of the focus S 
(in the general sense) of an optical system (Figure 3) , the latter being 
characterized by its exit diaphragm (real or virtual) . Point S is the parqxial 

image of a point — and monochromatic source S q . 

Let be the plane parallel to n^, with origin S and with position vector P. 
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Figure 3» Fraunhofer’s Diffraction. 
The diffracting plane ll 0 is the exit 
pupil of an optical instrument. The 
observation plane is normal to the 
axis at S, which is the paraxial image 
of a monochromatic point source S Q * 


To study the amplitude g(P) 
diffracted on the plane fFby pupil 11^, 

we use form (1. 6) of Huygens' 
principle which, as we have just seen, 
relates to the Fresnel diffraction at 
a finite distance. However, the 
intervention of the optical system 
modifies this expression. 

In effect it is necessary to 
introduce in (I. 6) the distribution of 
amplitudes h^M^ on plane n^; the 

distribution is given by the amplitude 
which would be produced on by the 


perfect spherical wave S which con- 
verges in S, in the absence of diffraction (since the diffraction intervenes only 
after passage through screen n^) ; whence, if r^ is the radius of S, 


( 1 . 12 ) 



By means of a reasoning similar to that of Paragraph I. 2 we obtain, in 
vectorial notation, 


(1.13) 


r = M S =-d + (m 2 / 2d 
o O \ 0 


Hence , after taking the independent term of outside the integral and 

taking (1.12) and (1. 13) into account, formula (1.6) becomes, save for the con- 
stant factor -i/Xd, 


(1.14) 


g(P) = e 


ikP 2 / 2d 



In this calculation we have assumed that the transparency of pupil is 

uniform. If this transparency is not uniform but represented by the complex 
function f^M^, (1.12) becomes 


(1.15) 
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This situation occurs especially in the filtration of spatial frequencies 
where a photographic plate of transparency t = f^M ^ is placed against il^. 


Under these conditions we obtain, instead of (1. 14) , 


(1.16) 


g(P) 


ikP 2 /2d 

e 



M dM 


which, save for the factor 


(1.17) 


ikP 2 / 2d 
e 


is the Fourier transform of f 



To obtain (1. 16) we have started out from formulas relating to the 
Fresnel phenomena, but the introduction of the term (1. 12) , due to the sphericity 
of wave 2, gives a Fourier’s transform as a final result save for factor (1. 17) . 


a. The Quadratic Phase Factor 


As the case of the Fresnel diffraction (Paragraph 1.3) the factor 
(1. 17) is frequently neglected — because it disappears during the recording of 
the phenomenon (1. 16) by means of a cell or a photographic emulsion. In 
coherent optics it is no longer permitted to ignore it; in effect: 

1) In Fourier-transform holography (1.17) is recorded, and 
it changes the properties of the hologram. We shall return to this point in 
Chapter IV; 

2) In the filtering of spatial frequencies it is necessary to 
compensate (1. 17) in order to obtain linear relations invariant by translation. 


b. Reciprocity of the Fourier Transformation 


The presence of the quadratic phase factor (1. 17) (which 
represents a spherical wave) means that the Fourier transform of function 
f^M \ is not formed on plane JF but on a sphere of radius d, with its center in O. 


Factor (1. 17) will often be compensated below (Chapter III) ; the distri- 
bution (1. 16) of the amplitudes on plane JF then becomes the Fourier integral 
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We know that this transformation is reciprocal and that we have 


(1.19) 



J F (P)exp(f 


P.M 


dP . 


c. Reduced Coordinates 


Let us write ( Figure 4) for coordinates x and y of point P of 
the Fourier plane 3\ 

(1.20) u = x/\d, y = y/\d; 


where u and v are the reduced coordinates of P which define the vector £2 (u, v) 
of the Fourier space. 


Hence relations (1.18) and (1.19) become, save for a constant factor, 


( 1 . 21 ) 


( 1 . 22 ) 



which corresponds to the usual notation of the Fourier transformation. 


d. Notations for the Fourier Transformation 

We shall generally use for a pair of Fourier transforms the 
same pair of letters, one of them being a lower case letter and the other one 
a capital letter. 

We shall also write, for the Fourier transformation (1. 21) , 

(1.23) f(M o ) F<T » F(S2) 


( F. T. = Fourier transformation) 



Figure 4. The Reduced Coordinates. 
If x and y are the coordinates in the 
Fourier plane, the reduced coordin- 
ates u ss x/Xd, v = y/\d define the 
vector fi(u,v) of the Fourier space. 


and for the inverse Fourier transformation 

F. T. -1 f^M ^ . 

(1.24) F(J2) ► \ 0 / 

Notes . 

a) The time factor chosen is e 1C0 ^; hence the phase lags are positive. 

b) In line with the usage adopted in optics we shall frequently use the 
term "amplitude" for "complex amplitude. " 


5. Conclusion 

A large number of diffraction problems may be solved by the 
approximate formulas (1.9) and (1.16) relating to the Fresnel and Fraunhofer 
diffractions. 

When the diffraction angles are large, formula (1.2) gives satisfactory 
results within the limits of the scalar approximation (study of the aberrations 
of the holograms) . 

Finally it should be remembered — and this is often useful in certain 
reasonings (for example, Paragraph II. 7) — that Huygens' principle may be 
strictly expressed, with the aid of a Green function, by a convolution. 


II. HOLOGRAMS 


1. Introduction 

Holography has introduced considerable progress into coherent 
optics. While originally the purpose of the holograms was to reconstitute two- 
or three-dimensional images from a photographic recording made without 
lenses, their present field of application is much more widespread. 
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Hence the term "lens-less photography" gives a rather restrictive 
meaning to this method, and it is more accurate to use the expression 
"recording of a coherent optical wave in phase and in amplitude. " 

In effect it is here that the interest of holography resides. By means 
of a coherent background it is possible to record any diffraction phenomenon 
on a photosensitive material. This phenomenon may then be restituted at will 
by coherent illumination of the photographic emulsion. In particular, if the 
diffracted wave originates from a three-dimensional object, the image of this 
object in relief is reconstituted, but in the case of a randomly chosen wave its 
manipulation is the equivalent of that of a photographic plate. 

Another remarkable property of holograms is that it is possible to 
record, successively, several waves on the same emulsion. In that case the 
process taking place during the restitution is an addition not of the intensity but 
of amplitude and phase. 

The principle of holograms is now sufficiently well known to enable us to 
take up its study in the form of a general method which is frequently fruitful in 
research. 

We shall stress more specifically the new aspects obtained herein: 

a) Properties of holographic image production by a simple geometric 
method; 

b) Diffuse illumination interferometry; 

c) Duplication of holograms; 

d) Applications to the filtering of spatial frequencies (Chapter IV); 

e) Experimental results in three-dimensional holography. 


2. Physical Aspect of Holography 

The idea which led to holography [2-4] was based on the Huygens- 
Fresnel principle and on the properties of luminous interferences. 

Let us take at random a diffracting system which emits a coherent 
electromagnetic wave. This system is represented schematically in Figure 5 
by an object O illuminated by a laser <£. 
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Figure 5. Huygens’ Principle. The 
radiation received by point P may be 
considered as originating either from 
object 0 illuminated by laser £ , or 
from a random virtual plane H situated 
between 0 and P. 


Let us take, in addition, a point P which receives the radiation emitted 
by O, and a virtual plane H situated between O and P. The wave emitted by 
object O passes through plane H before arriving at P. Thus it can be said that 
the radiation received by P originates from the object, that it can likewise be 
said that it originates from the virtual plane H. 

Let us now consider the following two situations: 

a) The object illuminates point P and gives a distribution of complex 
amplitudes A(x,y,t) on plane H(x,y); 

b) The object is absent, but a synthetic source has been reconstituted 
on H; the vibrational state of this source is A(x,y,t) . 

According to the Huygens -Fresnel principle the effect at P is the same 
for both situations. 

A hologram is precisely a photographic recording which, when illuminated 
in a suitable manner, restitutes the wave A(x,y,t). 

However, we know that it is not sufficient to place a photographic emul- 
sion in plane H in order to obtain such a recording. In effect the progressive 
monochromatic wave emitted by O will give, on plane H, the field distribution 

(II. 1) A(x,y)cos[wt - <p(x,y) ] , 

but the blackening produced on the plate will be a function only of the received 
energy, proportional to the mean, in time, of the square of (II. 1) , i. e. , to 

(II. 2) [A(x,y)] 2 

and the phase information will be lost. 
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To preserve both the phase and the amplitude information (II. 1) is 
transformed into a stationary wave by the addition of a coherent background. 

In this way we obtain a network of interference fringes whose recording con- 
stitutes the hologram. The local variation of the contrast of the fringes and of 
the fringe interval in the inter ferogr am then reflects the local variation of the 
amplitude and the phase of the wave (II. 1) . 

Experimental Arrangement. The coherent background is added simply 
by means of the assembly whose schematic outline is shown in Figure 6. 



Figure 6. Recording of a Holo- 
gram: Principle of the Assembly. 

A laser beam £is separated into 
two waves (1) and (2), by the 
separating L. Wave (1) illumin- 
ates object 0. The reference wave 
(2) directly illuminates the photo- 
graphic emulsion H; it constitutes 
the coherent background which 
interferes with wave (3) diffracted 
by the object. The plate records 
the total illumination due to waves 
(2) and (3). 


The coherent background is also 
called carrier or reference wave. 

Let us examine this process in 
greater detail, as well as the manner in 
which hologram H restitutes the recorded 
wave. 


3. Principle of the Holograms 

The experimental arrange- 
ment of Figure 6 is reproduced in a 
schematic manner in Figure 7. The 
same laser — not represented — 
illuminates the object and gives the 
reference wave. The object may be 
illuminated in any manner; for example 
one may interpose a polished glass 
between the laser and the object. 

The reference wave is generally 
a parallel, divergent or convergent, 
wave; hence it may be represented by a 
point source S . 


a. Recording of the 
Holograms 

We use the subscript 
o for the object and p for the carrier. 
The coordinates on H are x and y. 
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Figure 7. Recording of the Hologram 
(Shown Schematically). A laser 
illuminates the object in any given 
manner (if necessary through a polished 
glass). A part of the laser beam serves 
for forming the point reference source 
Sp . The hologram is recorded on H. 


The object diffracts on H — neglecting the time factor e — the com- 
plex amplitude 

(II. 3) A Q (x,y) = a o (x,y)exp 

This distribution is unknown a priori, but is perfectly well defined because the 
illumination is coherent. 


i(p o (x,y) . 


The carrier gives, on H, the complex amplitude 


(II. 4) 


A (x,y) = a (x,y)exp 

r r 


i<P p (x,y) 


Plate H records the total illumination due to (II. 3) and (II. 4) , or 
(II. 5) 


I = 

A + A 

2 = 

A 

2 + 

A 


o p 


P 


o 


2 + A A* + A " A 


o p 


o p 


The amplitude transparency t of the plate on which the image is produced is 
proportional, if y is the slope of the characteristic curve of the emulsion, to 
(Appendix 2) 


(II. 6) 



To free ourselves from the influence of the properties of the emulsion and to 
make the recording linear we bring about the following relationship (which can 
be easily done experimentally) 


(II. 7) 


A » 

A 

P 

o 


Then we can undertake an expansion of 
in mind, 


(II. 6); we obtain, keeping (II. 5) 
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(II. 8) 


2 


A* A 
y op 


f 


t - ia r Y n — , — 

p 


A 

y | o 

2 A 


>|c 

A A 

y o' p 


2 

A 

2 2 

A 


P 1 


P 


The reference wave is generally plane or spherical; hence, with the approxima- 
tions usually made in diffraction, it produces an illumination I A I 2 which is 
constant on H and which will be standardized as 1: ^ 


(II. 9) 


2 = a 2 = 1 


Condition (II. 7) makes it possible to neglect the term (y/2) 


Hence (II. 8) may be written, except for the constant 
a uniform absorption of the holograms r 


-y 


A /A 

o p 


which represents 


(II. 10) 


t = i -Ia a* -|a*a 

2 o p 2 o p 


In this way we obtain an expression which does indeed contain the ampli- 
tude and the phase of the diffracted wave, because (II. 10) can also be written as 


t = 1 - ya o (x,y)CQS 


cp (x,y) - xp 
o pj 


In this form the interferential aspect can be readily seen. It may be noted that 
(II. 10) is real and positive, and the complex recorded wave (II. 3) appears 
only by diffraction during the restitution. 


b. Reconstruction 


Hologram (II. 10) is illuminated by a reconstruction wave A^ 

which is identical to the reference wave (II. 4) which was used during the record- 
ing; this is not indispensable, as we shall see later, but the demonstration of 
this process is simplified. Thus on the minus face of the hologram (Figure 8) 
the amplitude is 

(II. 11) A - A . 

r p 

On the plus face of H the transmitted amplitude is A. = tA^, or 


( 11 . 12 ) 


A. = A 
l r 


4a a 


A 

o p r 


- I A* A A 
2 o p r 
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Figure 8. Reconstruction of the 
Hologram. The reconstruction 
source $ r illuminates the hologram 
H. The latter transmits three waves: 
1) the slightly altered reconstructed 
wave (cross-hatched); 2) two waves 
originating from two reconstructed 
images 0] and O 2 • The two images 

may be real or virtual, or partly vi rtual 
and partly real, depending on the 
position of S r at reconstruction. If S r 

is identical with the reference source 
Sp (Figure 7), one of the images is 

identical to the recorded object and 
is situated in the same region; 


This is the general expression of the amplitude transmitted by the hologram. 
Keeping in mind ( II. 9) and (II. 11) , it becomes 

(11.13) A. = A - ^ A -Za 2 A* . 

1 ' i r 2 0 2 r o 

Thus the hologram reconstructs three waves. 


1) The first term A^ of (II. 13) corresponds to the directly transmitted 

reconstructed wave. In reality this wave is slightly perturbed by the term 
-(y/2) IA q /A^ 2 which was neglected when going from (II. 8) to (II. 10) . 


2) The second term represents, save for the factor -"y/2, the ampli- 
tude A^ diffracted by object O on H during the recording; hence this wave is 
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wholly reconstituted in phase and in amplitude. An observer situated to the 
right of H will, in particular, see a three-dimensional image having all the 
characteristics of object O: relief, parallax and depth of field ( Figures 17 
through 26) . If object O is a wave surface given by any given coherent system, 
this wave surface will have been reconstituted. 

3) The third term, -(y/2) A 2 A* , represents a second reconstituted 

r o 

image 0 2 which will be called conjugated image, as distinct from image Oj which 
is called direct image. The properties of the conjugate image will be completely 
defined in the study of a point object and source presented in Paragraph II. 5. 

Let us note, for the time being, that this image may be virtual, real or partly 
virtual and partly real, depending on the operating conditions. 


4, Notes on the Holograms 

Several conclusions may be drawn from the reasonings just made. 

a) They of the emulsion (expression 11.13) intervenes only in 
the form of a constant factor in the amplitude of the waves diffracted by the 
hologram. The sign of y , which is + for a negative emulsion and - for an 
inversible emulsion, fixes the phase of the diffracted waves to within i r. Hence 
a copy of a hologram, which reverses the sign of y , has the same properties as 
the original hologram. This fact, in particular, makes it possible to use the 
method of hologram reproduction proposed in Paragraph II. 8. 

b) Contrary to general photography which consists of a point-by- 
point correspondence, a hologram results from the point- (whole) object 
plane correspondence of the hologram. Hence it has been proposed that any 
given portion of the hologram is sufficient to reconstruct the whole object. 

While this is true for diffusing objects which are then rendered with a decreased 
resolution and at a given angle, this is no longer true for nondiffusing objects. 

In particular , in Fourier holography (Filtering of the Spatial Frequencies, 
Chapter IV) a given portion of the hologram permits the reconstruction of only 
some spatial frequencies of the object. For a periodic object it is even possible 
that nothing is reconstructed from the whole. 

In interferential holography a portion of the hologram permits the obtain- 
ment of only a portion of the interferogram of the same order of magnitude. 

This has led us to develop an interferential method in diffuse light. 

c) The angular separation between the three transmitted beams 
is due to the use of an "inclined" carrier. In first approximation this means 
simply that the cone, defined by the carrier resting on the contour of the hologram, 
need not contain the object. Thus a separation can take place even when the carrier 
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is not inclined but is normal to the plate. The geometric study contained in 
Paragraph II. 5 makes it possible to treat these field problems in an easy 
manner. 


5. Geometric Study of the Holograms 

Now that we are familiar with the principle of holography, we can 
make a purely geometric study of the formation of the images by the holograms. 
The latter, in effect, have numerous analogies with lenses. 

The formulas proposed below seem to be easier to use than those given 
elsewhere in the literature [9,10 ] , In effect the axes of the different beams 
involved in the holographic technique are different ones; thus one has the 
advantage of using polar coordinates which simplify the results and their 
interpretation, 


We shall conduct our reasoning, as in classical geometric optics, on a 
point object M q which in this case gives two images M t and M 2 . The reference 

wave (carrier) is emitted from a point source and the restitution wave from 

a point S . Points M , S , and S have arbitrary positions; they may be real 
r opr 

or virtual, at a finite distance or at infinity. 


We shall examine successively the formation of the images in the 
meridian plane and outside the meridian plane. 

a. Formulas in the Meridian Plane 

(1) Notations (Figure 9), Points M^, S , and are defined by their 

"polar" coordinates (p,0 ) q The radius vectors p have their origin in any 

given point C of hologram H. The latter is situated in plane z = 0. 

(2) Sign Conventions . The p's are negative for z < 0 and 
positive for z > 0. The 0's are marked off as on the figure, in space z > 0. 

They are positive in the counterclockwise direction. The origin of the angles is 
axis Cz. 1 0 I is always < tt / 2 . 

(3) Establishment of the Conjugation Formulas . On recording, 

points M and S diffract amplitudes A and A on H. On restitution the source 
o p op 

S diffracts amplitude A on H. 
r r 
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P<0 P< o 



Figure 9. Geometric Study of Holograms; Notations, 
The object M Q , the reference source S p and the 

restitution source S r may be real (p > 0), virtual 
( p > 0) at a finite distance or at infinity. Arrow f 
indicates the positive direction for 0, 


It is sufficient to substitute these values into the general expression of 
the two restituted waves (expression II. 12) : 

(11.14) direct wave A i = A A' 1 ' A 

o p r 


conjugate wave A 2 = A^A^A^ 

in order to obtain the properties of the two images. 

Let us determine A , A , and A by doing our reasoning for the point 
o p r 

object M q (Figure 10) . Let C be the origin of the phases. Point M q radiates 
to point P, of ordinate y, on H, the amplitude 

(II. 15) A q = exp ik/r Q - p ^ . 

In triangle M^CP we have the following relation, keeping in mind the 

signs of 0 , p , and r : 
o o o 

r 2 = p 2 + y 2 - 2p y sin 9 , 

o o J o o 


whence if 



» 


lyl , and when 0 q is small, 
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H 


(II. 16) r o -p o -(5rV%> o ) 


and (11.15) becomes 


(II. 17) A o = exp ik/- - y0 \ . 


Likewise the carrier gives on H 


A = exp ik (— - yO 

p (2p p y p 


A_ = exp ik(~ - y« 


Figure 10. 

(11.18) 

and the restitution wave 
(H. 19) 


In these expressions coefficient 1/p of y 2 /2 is the curvature of the wave and 9 
is the angular coordinate of the center of curvature. 

Substituting expressions (11.17) to (11.19) into (11.14), the two restituted 
waves may be written as follows: 


Direct wave: 

y 2 

(II. 20) Aj = exptik 
Conjugate wave: 


r 1 i i 


2 \p p p 

\ o p r 


+ — - y 0 -0+0 


P r 


(11.21) A 2 = exp{ik 


y / 1 i 


2 l P n P v 

\ o p r 


+ _ + _ - y -0 + 0 + 0 


opr 


These waves converge to (or diverge from) their centers of curvature M t and M 2 
which therefore are the restituted images. Coordinates p and 9 of points Mi 
and M 2 are given by 


( 11 . 22 ) 


JL-JL _L _L J_ = i_+_l+i_ 

Pi P P P ’ P 2 P P P ’ 
o P p r o p r 
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(II. 23) 


01 = e - e + e , e 2 = -e + e + e . 

opr 4 opr 


direct image Mj conjugate image M 2 . 

These expressions are conjugation formulas analogous to those of 

lenses. They make it possible to determine the position of the direct image 

Mi and that of the conjugate image M 2 defined by coordinates p it 0! and p 2 , 0 2 , 

as a function of the coordinates of the object M , of carrier S , and of the 

restitution source S . 0 f’ 

r 


If pi is positive, Mi is real; if p\ is negative, Mj is virtual and the 
same applies to M 2 . All combinations are possible for the two images, as we 
have already indicated above. 

Furthermore we obtain formulas for the [angular] enlargement and 
[transverse] magnification of images Mj and M 2 . 

(4) Angular Enlargement . To define this enlargement we 
must specify the position of the eye. We shall assume — and this is logical — 
that the eye is situated in the plane of the hologram. By differentiating (11.23) 
we obtain 


(11.24) g! = d0j/d0 o = +1 , g 2 = d02/d0 o = -1 . 

(5) Transverse Magnification . It is given by 

0i,2 = d(p0)i )2 /d^p o ,0^ 


with p constant, or, keeping (11. 24) in mind, 


(11.25) 


(11.26) 


01 “ Pl/P Q > 02 “ P Q 


(6) Axial Magnification . 


dpi 

ffi = hT 1 = i^ 




From (11.26) it is inferred that, for a three-dimensional object, the direct 
image possesses a normal relief and the conjugate image an inverse relief. 
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(7) Change of Wavelength . If the recording of the hologram is 
carried out with wavelength X and the restitution with wavelength X’ , it is suf- 
ficient to substitute k f = 2i r/X' instead of k = 2ir/X into (II. 19) in order to obtain 
the conjugation formulas: 


(II. 27) M t 


_L = 171 
p, = x\p 0 



o 




Pi 





+ e 

p 


+ e 

r 


V v 

The enlargements and magnifications are simply multiplied by X'/X. 


b. Extra-Meridian Conjugation Formulas 


(Figure 11) . 


Points M , S , and S are now defined by coordinates p , 6, an ^ <P 
or p 



By means of a reasoning similar 
to that used in Paragraph II. 5. a we 
obtain conjugation formulas which are 
still the same as (11.22) and (11.23) for 
p and 0. To these are added the for- 
mulas containing cp: 


(11.28) 


(p i = cp - <p + cp , 
^ 1 o p r 


cp o = -cp + 0 + cp . 

^ 1 o p r 


Figure 11. Configuration for the 

Extra-Meridian Conjugation Formulas. For the magnifications nothing 

has changed; as for the enlargements, 
they must be completed by the formulas 
containing cp : 


d (p i/ dcp = 1, 

(11.29) 

d<p -J dip o = -1. 
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Invariance of the Formulas. The position of origin C on the hologram 
has been chosen at random for establishing the conjugation formulas; hence 
the latter are invariant when C is displaced. 


Likewise the origin axis of the 0's may be chosen at random. In par- 
ticular, if instead of the z-axis we take the axis S C of the restitution wave 

r 

for origin (Figure 9) , formulas (II. 23) become even more simple: 


(11.30) 


0i = 9 - 9 , 
o p 




c. The Aberration of Holograms 


For the establishment of the conjugation formulas we have 
neglected, in expression (II. 16) , all the terms of an order higher than 2. 
This is comparable to the Gaussian approximation in classical optics. The 
higher -order terms represent the aberrations studied by various authors 
[9,11], In this way one obtains the usual aberrations, spherical aberration 
coma, distortion, astigmatism, and "chromatic aberration" in the case of 
a change of wavelength on reconstruction. 


We note that it is possible to obtain the conditions of rigorous stigma- 
tism for the two images without long mathematical expansions. To this end 
we start from the rigorous expressions of the amplitude of the reconstructed 
waves which are, according to (II. 14) , 


(11.31) 

(11.32) 


Aj = A A' 1 ' A = exp 
opr 


A 0 = A ' A A - exp 
6 opr 


i 1<P - <P +(p 

\ o P r 


il-cp + 0 +• w I 

\ o p r 


(1) Direct Nonaberration Image . Let us take a point object 
M q and its direct image M* (Figure 12) . Stigmatism prevails if, regardless 

of the position of point P on the hologram, the optical path l M^PMj J is con- 
stant, or, according to (11.31) , ' / 


(11.33) 


CD -CD + CD 

o p r 



Constant 


where cp , (p 
r p 


constant, which entails, save for a constant factor , 


(II. 34) 


A = A , • 

r p 
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(11.35) Ax = A . 

Hence rigorous stigmatism prevails if the recon- 
structed wave is identical to the reference wave; 
the latter can be a deformed wave taken at random. 
Then the image is identical to the object and is 
situated in the same region. 

In the special case of a spherical reference 
wave originating from a point source S^, the condi- 
tion of stigmatism is that the point reconstruction 

source S must be identical with S . 
r p 

(2) Conjugate Nonaberration Image. In 
the same manner we find that, if 

(11.36) A = A* , 
r p 

the conjugate image is stigmatic. Moreover, (11.32) 
is then written as 

(11.37) A 2 = A* . 

These formulas signify (Figure 13) that the nonaberration conjugate image 
0 2 , the corresponding reconstructed source S ^ and the same elements 

relating to the direct image are symmetrical with respect to H. 

Thus if S , is real, S „ is virtual, and vice versa, and the same 
rl r2 

applies to the images. 

d. Conclusions Regarding the Production of Holographic Images 

The conjugation formulas immediately give a number of general 
rules, the most important of which are cited below: 

(1) Angular Dimensions of the Images . According to (11.24), 
regardless of the operating conditions, the two images are seen from a point 

C of the plane of the hologram at an angle that is always constant, which is 
that at which the object is seen from the same point. 

(2) Position of the Images in Space. According to (11.30), if 
we take axis CS^ of the reconstructed wave as the origin of the angles, the 



Figure 12. Rigorous 
Stigmatism. M Q is the 

object; Mj the image. 

The condition of rigorous 
stigmatism is, according 
to Fermat’s principle: 
M 0 PMi = constant. 
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Figure 13.- (a) Conjugate waves. Two conjugate luminous distributions 
A and A*, on H, correspond to two waves Sand S* which are symmetrical 
to H because to go from Sto 2 * one simply changes the sign of the path 
difference A. (b) The nonaberrant positions. (1) Direct image M] : the 
reconstruction source S r i must be identical with reference' source Spj 
then image M] coincides with object M 0 . (2) Conjugate image M 2 : 
restitution source $ r 2 is symmetrical to Sp with respect to H. The con- 
jugate image is then symmetrical to the object. 


directions of the two images Mj and M 2 become 
(11.38) &i = - 0 , 0 2 = -0J. 


Consequently these two directions are symmetrical relative to the 
reconstructed axis (Figure 14) and are invariably linked to the latter by angle 
0 q - 0^ between the carrier and the object, defined during the recording. 

Hence, 

a) If the reconstruction source turns by an angle a , the two images 
turn by the same angle; 


b) 

space. 

write p = 
P 

(11.39) 


If the hologram turns by an angle a , the images remain fixed in 


(3) Plane Reference and Reconstructed Waves . In (11.22) we 

p = 00 , whence 
r 

Pi = -P 2 = P Q • 
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Figure 14. Position of the 
images in Space. 


The two images are linked in an invariable man- 
ner to the restitution axis S^C by the angle 

between carrier and object, 9 - 9 , defined at 

0 p 

the time of the recording. 

If the reference wave and the reconstructed 

waves are plane, the direct image O t is identical 

to the object, and the conjugate image is deduced 

from O t by symmetry relative to a plane whose 

trace is the CN normal to S C. 

r 


Hence, according to (11.38) , the two images are symmetrical with 
respect to a plane whose trace on Figure 14 is the normal CN to the recon- 
structed axis. 


(4) Limitations of the Formulas (Bragg's Effect) . One of 
the limitations of the use of the conjugation formulas is due to the aberrations. 
Another limitation is of a physical nature. In effect the emulsion is not 
infinitely thin as has been assumed until now. Its thickness is finite, of the 
order of 15 m in the case of the Kodak 649 F plates. 

Thus the plate records the state of interferences not on a plane but in 
the volume occupied by the gelatin. We then obtain silver layers which con- 
stitute a three-dimensional network. 


During the reconstruction this network presents an effect which is 
more directional when the fringes are closer to one another'and more inclined 
in the emulsion. 

Hence if during the recording - ,0 p and(u0 Q + 0 \/ 2 are small, the 

events taking place conform to the formulas, but when these angles are large, 
the diffracted amplitude rapidly decreases when 9 ^ deviates from 0 p and the 

image disappears. 


Example . For the hologram of Figure 19, we had 9 ^ - 0^ - 35°; if 

the reconstructed wave 6 deviates by ±10° from angle 9 of the carrier, the 

r p 

diffracted intensity drops to 10 percent of the maximum value obtained for 
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Figure 15. Assembly for the Recording of Holograms, with Two Illuminating 
Beams. The beam originating from a laserJSis divided into three parts by 
plates L ] and L_2- Two of the beams serve for illuminating the object 0, and 
the third (cross-hatched) constitutes the reference wave. All the beams are 
dilated by the lenses of microscope M in such a way as to cover plate H and 
the object. The reference source is then the focus Sp of objective Mp . 


Furthermore, for 9 ^ ~ 0 , the conjugate image is extremely weak; 
to observe it, it is necessary to give 6 a value close to 0 . 


6. Experimental Study of the Holograms 

This paragraph deals with three-dimensional holograms. Inter- 
ferential holography is studied in Paragraph II. 7. The two-dimensional 
Fourier holograms will be examined in Chapter IV. 
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a. The Assemblies 


(1) Recording . The assemblies using parallel light exhibit a 
disadvantage of a practical nature: the beams are produced by collimators 
which are bulky and, in addition, the two-dimensional waves cannot be filtered 
from their heterogeneities. Hence we have used the assembly using divergent 
light, as shown in Figure 15. 



Figure 16. Reconstruction 
Assembly. The hologram is 
illuminated by a monochro- 
matic point source S r pro- 
duced by laser £. 


Similar assemblies have been used with one, 
two, or three illuminating beams. 


(2) Reconstruction . The recon- 
struction takes place visually or photograph- 
ically as in Figure 16; point is in a posi- 
tion close to that of point of Figure 15, 


in order to prevent aberrations. Figure 17 
shows the whole restitution assembly: the 
image which appears in the hologram is that 
of the radar, shown also in Figures 21 
through 23. 



Figure 17. The Reconstruction Assembly. The image which appears in the 
hologram is that of the radar shown also in Figures 21 through 23. 
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Figure 18. Effect of Field Depth. The camera is focused on the foreground 
of the reconstructed image. The church in the background is blurred 
(Figure 19). 



Figure 19. Effect of Field Depth. The camera is focused on the stained- 
glass window and the shadow of the statue; only this part of the church is 
in sharp focus. The statue in the foreground is completely blurred 
(Figure 20). 
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Figure 20. Effect of Field Depth. The lens of the camera is stopped down 
to a diameter which is smaller than that of the eye; the whole scene is in 
sharp focus, but we note, in an accentuated manner, the granulation which 
is characteristic of diffusing objects illuminated by a laser and observed 
visually. 


For the reconstruction of Figures 18 and 19, the diaphragm of the lens 
was more open (30 mm) , the granulation disappeared, and so did the depth of 
field. 

b. Results 

Figures 18 through 26 show the reconstructions, obtained from 
holograms, of a number of three-dimensional objects. Note, in particular, .the 
parallax effect and the effects of the depth of field. 

c. Experimental Conditions 

The principal factors which influence the quality of holograms 
are the resolving power of the plates, the coherence length of the laser, the 
ratio between the amplitudes of the diffracted light and the carrier, and the 
mechanical stability of the recording arrangement. 
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(1) Resolving Power of the Plates . The photographic plate 
must record the microfringes of interferences which are formed in its plane 
( Figure 27) . 



Figure 21. Effect of the Depth of Field. Focusing on the ‘horn’ of the 
radar, on reconstruction (Figure 22). 
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Figure 22. Depth of Field Effect. Focusing on the parabola of the radar 
(Figure 23). 




Figure 23. Depth of Field Effect. Focusing on the shadow of the radar 
antenna. 


An exact study of the conditions under which the fringes are recorded 
requires the knowledge of the modulation transfer curve of the emulsion. This 
is generally not too well known. However, it is possible to do an approximate 
reasoning by means of which the orders of magnitudes are fixed. 


Let us consider two plane waves and 2^ originating from the object 

and from the reference wave, respectively (Figure 28) . Let us assume that 
and 2^ are inclined to hologram H by angles and -a . The hologram is 

then formed of parallel fringes whose interfringe 

(H* 40 ) Ay - A/2 sin o' 

must be resolved by the emulsion. 


Example , a = 30°, which is a common case; Ay = A. The resolving 
power of the plates should be close to the wavelength of the light. In practice 
for A = 6328 A (helium-neon laser) Kodak 649F plates are used, which 
resolve more than 2000 lines per mm. 

(2) The Coherence Length of the Laser. This limits the depth 
of the objects which may be recorded. In effect, let A and B be the extreme 
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points of the object (Figure 29), assumed to be illuminated according to the 
direction of the arrow. 

Then the path difference, relative to the carrier, cannot be equalized 
better than to within ± D. Hence it is necessary that the coherence length of 
the laser be greater than 2D. 

The so-called "monomodal" lasers actually possess several axial modes. 
The 3-mW helium-neon laser, which has served for the preparation of the 
holograms of Figures 17 through 26, had a cavity length of 1 = 75 cm, which 
corresponds to an intermode distance Av = c/21 = 200 MHz. In actual practice 
we have observed, as a result of beats, 3 modes, whence 2D = c/3Av = 50 cm. 

For long exposures, however, it is prudent to take into consideration 
the whole width of the Doppler ray, or 1000 MHz. 

Note that in Fourier holography, the path differences are small and no 
special precaution must be taken when using the laser. 

( 3) Amplitude of the Reference Wave . The photographic 
emulsions are generally not linear, which means, in the case of holography, 
that if the hologram receives ( Paragraph II. 3) an. illumination 

(11.41) 1= I A 2 + |A 2 + A A* +A* A 

o p o p o p 

its transmission factor in amplitude is not proportional to (II. 41) but to 

(11.42) t = f Y/2 . 

This expression may be expanded if 

(11.43) |A o /A p |«l; 

it becomes, when only the first terms are retained, 

A 2 A A* A* A 

o Y o P Y 0 P , 

A 2 A 2 “ 2 |A I 2 

P P I P | 

Physically this is the equivalent of choosing, on the blackening curve, an 
operating point determined by |A I 2 , and of modulating the curve by |A I 2 over 
a short length. P ° ' 
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Figure 24. Parallax Effect. This figure is the reconstruction of a set of 
dice and dominoes placed on a mirror. The images of the objects in the 
mirror are also reconstructed and can be differentiated from the object by 
their split appearance. In effect the mirror employed is silvered on its 
back, and furnishes two images. Focusing on the foreground (Figure 25). 



Figure 25. Parallax Effect. Enlarged detail of Figure 24; the lateral 
faces of the dice on the left are not visible (Figure 26). 




Figure 26. Parallax Effecto The camera was displaced to the right, the 
faces hidden in Figure 25 appear; the blurred dominoes in the background 
are also displaced. 



Figure 27, Photomicrogram of a Hologram. The interval between the 
interference fringes is of the order of 1 \x for this hologram. 
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Figure 28. Two plane 
waves Z Q and Sp inclined 
to plane H of the hologram 
by angles a and -a, produce 
on H interference fringes 
spaced Ay = A./2 sin a apart. 



Figure 29. Influence of the Coherence Length 
of the Laser. If object 0, illuminated by a 
laser in the direction shown by the arrow, is 
recorded by holography in plane H, the coher- 
ence length of the laser must be greater than 

2D. 


If condition (II. 43) is not satisfied, 
terms of a higher order of the type 


(II. 45) 



appear, which have two principal effects : 


a) They yield images of a higher order, like the ordinary networks. 
However, in view of the Bragg effect, these images disappear when the angle 
between the carrier and the diffracted light is large. 

b ) They change the distribution of the amplitudes in the two waves 
of the first order: the direct wave and the conjugate wave. This latter effect 
is not very pronounced in the three-dimensional holograms where the observa- 
tion is often visual. Nevertheless in spatial filtering where the hologram is 
made to fulfill a mathematical function (Chapter IV) , the result can be damaging. 
It is necessary in this case to respect condition (II. 43) . 


(4) Mechanical Stability of the Holographic Arrangement . The 
recording of a hologram consists, in fact, in arranging a photographic plate in 
the plane of an interferogram, in order to fix the latter. It is then necessary, 
especially if the exposure times are long, to take certain precautions in order 
to ensure the mechanical stability of the assembly. 


Let us use the example mentioned in Paragraph II. 6.c: the interfringe 
interval in the hologram was less than 1 ju: hence the production of a translation 
of the fringes of the order of 1 p, produced by the displacement of a single ele- 
ment of the recording device, is sufficient to make the hologram blurred. To be* 
sure, to obtain a high hologram quality it is necessary that the displacement or 
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the deformation of the fringes under the influence of mechanical vibrations or 
thermal effects be very much less than In, still referring to the above mentioned 
example. 

The installation of an antivibration support made of concrete, shown in 
Figure 42, has made it possible to avoid these effects. 


7. Interferometric Holography in Diffuse Light 
a. The Method 


The proposed method [12] makes it possible to prepare, by 
holography, an interferometer of the Mach-Zehnder type of a phase object, 
the latter being recorded through a diffusing screen. 

To this end the following two operations are carried out in either order: 

1) By means of the assembly of Figure 30 the hologram H of 

a polished glass D illuminated by a point source M is registered. The carrier 

is emitted from point S : 

P 

2) Without changing the assembly, the phase object O to be 
studied (Figure 31) is introduced in front of the polished object, and the image 
of the hologram of object O seen through the polished glass is prepared on the 
same emulsion. 

On reconstruction, the hologram on which two images have been pre- 
pared reconstitutes an inter ferogram, an example of which is given by Figure 32 
relative to a sheet of glass which exhibits "grains" and shows variations of thick- 
ness and of refractive index. 

It may appear surprising that this could be done through a diffusion 
screen; however, the explanation of the phenomena is simple. 

(1) Recording. In the first phase, point M gives on polished 
surface D the amplitude 

(11.46) a; 

The polished object transmits ( L is a linear operator) 

(11.47) L (a) . 

In effect the polished object does not destroy the coherence of the light, since 
its thickness is small compared with the coherence length of the laser. More- 
over Maxwell's equations are linear and the diffuser is passive, hence the action 
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Figure 30. Interferential Holo- 
graphy in Diffuse Light. First 
operation: The hologram of the 
polished glass D illuminated by 
the point source M is recorded 
on H. So is the reference source. 


Figure 31. Second operation: Without 
changing the assembly of Figure 30, the 
phase object 0 to be investigated is 
introduced between source M and the 
polished glass D, and the image of the 
hologram of 0 is prepared on the same 
emulsion. Thus the hologram of 0 is 
recorded through the polished glass 
(Figure 32). 


of the diffuser on the distribution (II. 46) is reflected by an unknown operator, 
because it depends on the nature of the diffuser; however, it is linear. Whence 
(11.47). 

Distribution (II. 47) gives, on the plane H of the emulsion (expression 
1. 9) the amplitude 

(II. 48) Aj = L( a) ★ d, 

where ★ signifies convolution and d is the normal derivative of Green's function 
relating to the diffraction (Fresnel's function in first approximation) . 

The carrier wave gives, on H, the amplitude 

(11.49) A . 

P 

The hologram records the illumination 

(II. 50) Ei = |A t + A p | 2 . 

In the second phase, the object O introduces the phase difference cp, and 
( II. 46) becomes 

(11.51) ae 1<P , 
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Figure 32. Interferential Holography in Diffuse Light. The hologram recorded 
according to Figures 30 and 31 reconstitutes an interferogram of the phase 
object 0 which has been recorded. In this case the object is a glass plate 
whose ‘grains’ and whose variations of refractive index and thickness are 
visible. The hologram obtained in this manner shows all the properties of the 
ordinary holograms in diffuse light, in particular, the interferogram is directly 
visible to the naked eye. 


Hence (II. 48) becomes 

(11.52) A 2 = L(a e 1C(> ) * d 


and the illumination (11. 50) becomes 
(11.53) E 2 = 


!A 2 + A 2 
2 P 


The total illumination, due to the two operations, is 
(II. 54) Ei + E 2 = I At | 2 + |A 2 1 2 + 2 I A 2 

[ P 

+ (Aj[ + A 2 ) A* + (A( + A 2 ) A^-. 
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(2) Reconstruction . Under the normal conditions of recording, 
the transparency t of the hologram (II. 54) is, according to expression (II. 10) 
in which y has been standardized at -2 

(II. 55) t = 1 + (Ai + A 2 ) A* + (A* + A* ) A . 

p p 

If (II. 55) is illuminated by the reconstructed wave A^ = A^, the second term, 
which reconstitutes the direct image, becomes 


(11.56) A t + A 2 

or, keeping in mind (II. 48) and (II. 52) , 

(II. 57) [ L(a) + L^a e 1 ^)] * d 

or, since operator L is linear, 


(11.58) 



This expression means that the reconstructed wave originates from the 
fictitious object 


(11.59) 


l[a (l + e 1(P ) 


which represents, physically, the projection, on the polished object, of the 
interferogram a (l + e^ ^of the dephasing object. 

b. Properties of Interferometric Holograms in Diffuse Light 


These are analogous to those of the usual holograms of diffusing 
objects. The grains or grooves of the emulsion are practically without effect 
on the quality of the interferogram. Moreover, to reconstruct the latter in its 
entirety it is sufficient to use any given part of the hologram. In actual practice 
1 cm 2 gives very convenient results. Finally the image may be observed by the 
eye directly. 


It may be noted that this is not the case in the methods using directed 
light [13, 14] : the reconstruction can be accomplished only by means of an 
optical system which entirely covers the surface of the hologram, since the 
area of the reconstructed image is directly proportional to that of the hologram 
used. 
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8- Reproduction of Holograms 

We shall conclude this chapter by discussing the important problem 
of the reproduction of holograms. This is beset with great practical difficulties 
[15] . In fact the interference fringes recorded in a hologram frequently have a 
spacing of the order of 1000 lines per millimeter; hence it is impossible to use 
optical systems for reproduction. Moreover, these fringes constitute a three- 
dimensional network (Paragraph II. 5. d) which rules out reproduction by contact 
printing. 

The process of coherent duplication proposed herein consists in using 
the beam directly transmitted during the reconstruction of a hologram H as the 
reference beam for a second hologram H', the object for H' consisting of the 
images of H. 



Figure 33. Coherent Copy of 
a Hologram H. When a plate 
H is used for the preparation 
of an image as shown on the 
figure, with S representing a 
coherent source of any given 
position, a hologram H’ is 
obtained which is an identical 
copy of the original, including 
the structure of the three- 
dimensional network recorded 
in the depth of the gelatin. 


Let us consider (Figure 33) a hologram 
H whose transparency in amplitude is, according 
to (II. 10), 

(11.60) t = 1 - -%r A A' - -%r A* A . 

2 o p 2 o p 

This hologram is reconstructed by illum- 
inating it with wave A^; it then transmits the 

amplitude distribution 

(11.61) A. = t A . 

l p 

Plate H' placed behind hologram H records the 
intensity 

(11.62) 1= |A.j 2 
and its amplitude transparency 


(11.63) 


t' = I 


-y/2 


becomes, using the same approximations as those 
which led to (II. 60) , 

(II. 64) t» a 1 + ^4 A A* + A* A . 

2 o p 2 o p 
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Except for some constants, this expression is the same as (II. 60) ; hence we 
have obtained a copy of hologram H. The quality of the reconstructed image is 
practically identical to that of the original, and the degradation due to the fact 
that the copy registers, in addition to its own background, that of the original 
as well, is not perceptible to the eye. 

It is shown, moreover, that plates H and H' need not be in contact, and 
that the copying wave may be chosen at random, as long as it is coherent and 
originating from a point source. 


9. Conclusions on Holography 

Holograms have two main functions: 

a) They permit the recording, on a photosensitive material, of 
both the amplitude and the phase of a coherent light wave; 

b) They permit the summate in phase and amplitude several 
coherent waves recorded successively. 

The first property makes the hologram a remarkable data storage sys- 
tem while the second represents a veritable data-pro cessing operation. 

In this chapter we have listed a number of other properties and applica- 
tions of holograms, but have reserved the applications to the filtering of spatial 
frequencies. The fact of being able to record an optical wave representing a 
complex function, e.g. , a Fourier transform, permits one, in effect, to obtain 
photographic transparencies which may act as complex filters in coherent 
optical filtering. 

Thus it is possible to develop certain applications of this method which 
had existed only potentially before the discovery of holograms. This will be 
discussed in Chapter IV. 


III. FILTERING OF SPATIAL FREQUENCIES 

1. The Possibilities Inherent in the Method 

Ever since its discovery in 1953 [1] , the filtering of the spatial 
frequencies in coherent light has been considered for a number of applications 
aimed, in particular, at the improvement of optical images [5-7, 16]. 
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This method has acquired an increased significance at the present time 
as a result of the new prospects which are opening up for this method in the 
optical treatment of data. This situation is due to several factors: 

a) Theoretical possibilities . A coherent optical system establishes 
between three functions of two variables — the signal f(x,y) , the percussional 
response h(x,y) of the system, and the response r(x,y) to the signal — a very 
general linear correspondence: a convolution relation. 

Now a number of mathematical operations are but special cases of 
convolution, notably the integration, differentiation, correlation, self- 
correlation or the scalar product of two functions ; 

b) Experimental possibilities . All these operations may be 
effectively carried out (as we shall see some examples of this in Chapter IV) 
not only for real and positive functions but also for complex functions. 

Experimentally, the operator "convolution by h(x,y) " is applied to 
f(x,y) in Fourier's space: this reduces simply to placing a filter into the 
spectral plane of the optical system. 

The gas laser and the holograms which permit the recording of com- 
plex functions greatly facilitate the realization of some of these filters or the 
material representation of complex functions. 

In this chapter we summarize the aspects of optical filtration which will 
be used later on, and give an example of application which may serve for the 
treatment of the photographic plates of bubble chambers. 


2. The Mathematical Bases of Optical Filtering 

Optical filtering in coherent light may be justified by purely 
mathematical considerations. 

In effect, the problem consists in finding the law of correspondence 
between a two-dimensional signal f and the (filtered) response r(M) given 

by an optical system defined by its percussional response h(M) . 

To this end it is sufficient to assume that 
/ 

a) the transformation is linear; 

^ b) the percussional response h(M) of the instrument is invariant 
by translation. 
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Condition (III. lb) means that if the signal is Dirac’s pulse , and 

if gives the response h( M) , then the pulse <5 , displaced by M* 

by translation, gives the response h(M - M*) » It is shown [17, 18] that condi- 
tions (III. 1) entail the law of correspondence. 

(III. 2) r(M) = f(M) * h(M) . 

This convolution product is written, after a Fourier transformation, as follows, 
with R(J2) , F( 12) and H(S2) being the Fourier transforms of r(M), f and 

h(M) , and S2(u, v) the general point in Fourier's space. 

(III. 3) R(fi) = F(S2) H<0) . 

From this it is deduced that the correspondence between f^M^ and r(M) is a 
linear filtering which follows the following scheme (® signifying multiplication) : 


(HI. 4) 


40 

Signal 


F.T. 

Spectrum 

analysis 



Filter 


The preceding reasoning does not make use of any physical mechanism. 

It is applicable to the case of both coherent and of incoherent illumination, 
provided the conditions (III. 1) are satisfied. However, in incoherent illumina- 
tion, the Fourier spectra F(I2) , H(S2) have no physical meaning. The case is 
different in the case of coherent illumination: in that case the scheme of 
( III. 4) may be realized by double refraction. 

3. Double Refraction 

a. Process 

The treatment (III. 4) may be carried out optically in a resti tutor 
by double diffraction [18]. The signal produces, by means of a first diffraction, 
a Fourier spectrum which, when filtered, gives (restitutes) the response by 
means of a second diffraction. 
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In reality the signal gives, in the spectral plane, the product of the 
Fourier spectrum and a spherical phase term ( Paragraph I. 4. a) . To obtain 
a linear filtering, i. e. , to accomplish the conditions (III. 1) , it is necessary 
to compensate this factor. This may be done by a convenient arrangement of the 
optical system [19] (also Paragraph IV. 3.b) . 

Then it is possible to neglect the phase term in the filtering process 
[Note: The compensation is effective only at the end of the treatment process, 
in the plane of the response (image) ; hence when a holographic recording is 
made on the level of the spectral plane ( Paragraph IV. 4. a) , it will be necessary 
to take into account these terms which are recorded in the form of holographic 
lenses, with their aberrations. ] 

This process consists of three phases which will be recapitulated below; 
the proofs will be given in Appendix 1. 

(1) The Restitutor . This consists of an object plane (signal) 

II , a spectral plane (F, and a response plane (image plane) II (Figure 34). 
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Figure 34. The Restitutor Operating by Double Diffraction. Object f M 0 is 
arranged in plane ll 0 ; it is illuminated by the convergent beam originating 
from lens L and gives by means of a first diffraction, in plane (F, the 
Fourier spectrum F(£2). The latter is filtered by a filter H(fi) situated in JF, 
and gives by means of a second diffraction, in plane II, the image of ll 0 

through lens L, the filtered response r( M) = f(M) ★ h(M). 
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The optical system consists, schematically speaking, of spectrum 
analyzer lens and of a restitution lens L. 

The signal to be filtered is generally recorded on a photographic 
emulsion whose transparency in amplitude t (Appendix 2) is proportional to 
signal f^M^ . 

(2) Spectrum Analysis . The signal f .M^ is placed in object 
plane having a general point M^x^, y^j . It is illuminated by a convergent 
beam originating from lens L^. 

The plane normal to the axis at the point of geometric convergence S 
of the beam is the spectral plane iF. The position vector in this plane is 
P(x t , yi) or, in reduced coordinates ( Paragraph I. 4. c) , 

(III. 5) S2(u, v) = P/\d, 

where d is the distance between planes and F, and X the wavelength. 

Under the effect of the convergent illumination, signal f ) diffracts 

a wave; the latter forms in plane F a luminous distribution which is the Fourier 
spectrum F(S2) of f^M^. 

(III. 6) f^I^ _£l!l_^F(fi) 

save for a phase factor. 

( 3) The Filtering . A filter is placed in plane F whose 
transparency in amplitude is H(fi) . Actually, in a limiting case, this filter 
may be a simple aperture. 

After passing through the filter the amplitude of wave F(S2) becomes 
(HI. 7) R(S2) = F(fi) H(fi) 

b. Restitution 

To ensure condition (III. 1) and to compensate the phase factor 
in (III. 6) , the restitution lens L should have its entry diaphragm in the spectral 
plane F. This is shown schematically on Figure 35 by a thin lens, placed in 
plane F. The transverse magnification of L will be standardized at-y = +1. 
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Figure 35„ Restitutor Operating by Double Diffraction. To ensure the 
conditions of linearity and invariance by translation, the diaphragm of 
restitution lens L must be placed in the spectral plane 5. 


After passing through L, the wave (III. 7) gives by means of a second 
diffraction the following response on image plane II having a general point 
M(x,y): 


(HI. 8) r(M) = f F(S2) H(fl) e 27ri R,M dO 

If. 

which is the inverse Fourier transform of the filtered spectrum: 

F. T, -1 

(IIL9) F(fi) H(fi) — > r(M) . 

Let h(M) be the percussional response of filter H(fi) : 

H(12) — - T ‘ . h(M); 

which represents, physically, the distribution of the amplitudes in the image of 
a point luminous object placed at the origin of plane 'll . 

Applying the convolution theorem to (III. 8) we then obtain 

(HI. 10) r(M) = f(M) * h(M) 

To sum up, we have: 


£ ( M o)' 

signal 

(III. 11) 


F.T. 


•F(fi) — R(fi) = F(fi) H(0). 


F.T* 


spectrum 


filtered spectrum 


H(Q)- 

filter 


F. T.' 


-l 


-r(M) = f(M) * h(M) 
response 


«.h(M) 
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4. Properties of Spatial Filtering 
a. The Spatial Frequencies 

The spectrum analysis correlates with each elementary sinusoidal 
component of the object f^M^ with point S2 of the spectral plane, and the 

amplitude of the light wave at this point is proportional to the amplitude to the 
elementary sinusoid. 

In effect, let us take such a sinusoidal component f^M^ of period p, 

and a unit vector U ( Figure 36) normal to the "lines" of the sinusoid. In a 
complex form we have 

(III. 12) f^M^ = a exp ^2 it i -£■ . M q 

We say that the spatial frequency of the object f^M^ is 1/ p or, to take 
the orientation of the object into account, U/ p. 




Figure 36. The Spatial 
Frequencies. The object f 
(M q ) = a exp |2 7ri (U/p).M^j , 

of spatial frequency 1/p, 
gives in spectral plane £F, a 
spectrum which, when f(M Q ) 

is unlimited, is a point P 
whose position in reduced 
coordinates is fl = U/p. If 
the object is not unlimited, 
this point is the center of a 
diffraction spot. 


In the spectral plane we obtain the 
Fourier transform of (III. 12) 


F (12) = a f exp 


2 7riM .( - - 0 
° \ P 


dM 


which is the Dirac distribution F (ft) 

= a 6 (q - — ^ . This represents a point 


with amplitude a and situated at 0 = 


_U 

P 


i. e. , at a distance from the center of the 
Fourier plane which is equal to the spatial 
frequency 1/p. 


Coming back to the simple coordi- 
nates in the spectral plane, with P= \d£2 
we have 


F(P) = a 6 
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in other words a point of the spectral plane fF situated at P = (Xd/p) U whose 
distance from center S is proportional to the spatial frequency and to the dis- 
tance d between planes II and 3\ 

o 

Thus the center of the spectral plane corresponds to the continuous 
component of the object and the edge at the high frequencies. 

NOTE. If, instead of considering the expression (III. 12) , which is the 
result of a decomposition of the object into complex sinusoidal components, 

we had considered the real sinusoid f (m Y= 1 + cos [ 2 7r — . M J , we would 

\ oj \ P oj 

have obtained the frequencies £2 - 0, corresponding to the continuous component , 
and £2 = ± U/p corresponding to the. decomposition of the cosine. 

b. The Cutoff Frequency 

If the filtering function H(£2) reduces to a simple aperture 
(the entry diaphragm of L, of transparency 1 and of radius R or R/Xd in reduced 
coordinates) , it only permits the passage of spatial frequencies such as I £2 1 
= 1/p < R/Xd. The spatial frequency 1/p = R/Xd is the cutoff frequency. 

The period of the corresponding "cutoff" sinusoid is seen from the 
center S of plane fF at an angle a. = p/d = X/R, which is the classical "resolving 
power" (in coherent light) due to the diffraction. All the spatial frequencies 
such as l/p < R/Xd are transmitted without alteration, since H(£2) =1 in the 
diaphragm; the instrument is a low-pass filter with a flat band. 

Below the attention is devoted to the data-treatment aspect and not to 
the power of resolution. Hence it will be assumed that all the useful frequencies 
are transmitted by the diaphragm. 

This means, in particular, that if H( £2) = 1 we have (expression III. 9) : 

r(M) - f(M) . 

The response is identical to the signal. 


5. An Example of Application 

To conclude this general part of the discussion it may be useful to 
give a simple example of optical filtering, which may find an application in the 
treatment of photographic plates of bubble chambers. 
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Figure 37. Extraction of Straight Lines by Optical Filtering. When 
the Fourier spectrum of the object (above) is filtered by a fine slit 
passing through the center'of the spectral plane, there remain in the 
image plane of the restitutor only the straight lines normal to the 
slit. By turning the slit, all the lines of the object are successively 
extracted. By contrast, if the filter had consisted of ao opaque line, 
the lines of the object normal to the filter would have successively 
disappeared ■ 

These plates comprise essentially parallel traces produced by the 
exciting particles. It may be advantageous to "cleanse" the plate of these 
traces, or on the contrary, to bring them to light. 

This can be done quite simply, because a long and fine straight line L 
gives, regardless of its position in the object plane, a spectrum concentrated on 
a line L' normal to L and centered on center S of the Fourier plane. 

To eliminate L and all straight lines parallel to L, it is sufficient to 
eliminate L' ; to ensure that only lines L will appear it is sufficient to allow 
only L’ to pass through. Figure 37 shows the result obtained when a slit is 
rotated in the spectral plane: straight lines L normal to the slit successively 
appear in the image. 
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IV. PATTERN RECOGNITION BY FILTERING OF SPATIAL FREQUENCIES 


1. Introduction 

Pattern recognition consists in detecting the presence and localizing 
the position of a signal in the presence of other signals or in the presence of 
noise. An example of such a process is the identification of characters or 
words in a text. 

It may appear that an optical self-correlation carried out on the signal 
would furnish a convenient means for accomplishing the identification. 

However, the methods using self-correlation lead to ambiguities, and 
it frequently happens that the object identified is different from that sought. 

To remedy these disadvantages a method has been proposed which is 
characterized by a low probability of false alarm [20, 21] . 

To this end the object to be treated (signal to be identified + noise) is 
subjected to a transformation which permits the survival of only the lines of 
discontinuity of the object, or more generally, the zones which exhibit a rapid 
variation of amplitude. 

The object transformed in this manner is then correlated with a signal 
model which had undergone the same transformation; the result of this is an 
identification with a high signal/ noise ratio. 

The transformation is obtained by the use of suitable differential oper- 
ators. The complete identification treatment and, in particular, the differen- 
tiation operations are easily carried out in coherent light on the basis of the 
properties of the Fourier transformation and of the filtering of spatial frequencies. 
It may be noted that this is a general method, and that it is just as easy to 
construct a filter leading to the identification of a complex object as a filter 
identifying a simple object. 

2. Principle of Pattern Recognition 

Let us consider a two-dimensional signal s in the presence of 

an additive noise b(M J. Using the notations of Chapter III, this corresponds 
to the object 



The problem consists in detecting the signal in the object. This may be 
done by filtering the spatial frequencies according to the scheme of Paragraph 
III. 2 as follows: (Note: When there is no ambiguity, the subscript "o" relating 
to the object will be dropped. ) 

/ 


(IV. 2) 


object 

spectrum 


f (M) = s(M) + b(M) 
F.T. 

F(fi) = S(f2) + B(f2) 
I 

— H(fl) filter 
I 


filtered spectrum R(£2) = S(fi) H(£2) + B(J2) H(ft) 

F. T." 1 

response r(M) = s(M)"fr h(M) +b(M)^ h(M) 


Response r(M) of the system consists of two terms. The first term is 
the convolution product r = s(M) ★ h(M) , which is the response in II due to 
signal s(M). 


The second convolution product r^ = b(M) ★ h(M) is the response due 
to the noise b(M) . 


In general, r g and r^ are no longer images but distributions of luminous 
amplitude. 

In actual practice, the identification of the signal s(M) consists in 
finding a filtering function H(fl) which would furnish in the image plane n 
(Figure 34) a distribution of illuminations E = |r (M) I 2 whose maximum E 
indicates, without ambiguity, 

a) The presence of the signal in the object; 

b) The position of the signal ( Figures 38 and 47) . 


This implies that the signal/ noise ratio must be sufficiently high. 

It is possible to determine from scheme (IV. 2) — without for the time 
being specifying H(S2) — some general properties of pattern recognition. 

a. The Signal/ Noise Ratio 

Let us assume that function H(ft) has been determined. H(fi) is 
the Fourier transform of the percussional response h(M) of the optical system 
f provided with filter H(fi) ] : 
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h(M) 


F, X -. H(fl) 

According to (IV. 2) the response in plane Ilis 

r(M) = s(M) ★ h(M) +b(M)^ h(M). 


Let |r I 2 and |r, I 2 be the maximum illuminations which 
s m b m 

responses. The signal will be detected if |r I 2 > |r, I 2 , 

. ./ . & .. | s |m I b lm 

if the signal/ noise ratio 


correspond to these 
or more precisely, 


(IV. 3) 



is sufficiently large so that the contrast c = ^ |r g 1^ - |r^ l^j/ ^Ir^ 1^ + |r g 
can be measured. Figures 55 and 56 give examples of such measurements. 



b. Automatic Positioning of the Response 


The properties of (III. 1) of linearity and invariance by transla- 
tion have two consequences. 

(1) Translation. The filter, designed for the recognition of 
signal s(M) , will also recognize the signal after a translation of Mj in its plane; 
the response in plane Ilis simply translated by the same quantity (Figure 
38). 

(2) Additivity 
give the response Sn r g ^M - M 

Property (2) is obvious; property (1) can be readily proved. 

Let us take a signal s(M) situated, for the sake of convenience, in the 
center of the object plane (Figure 38) . 

Let 6 be the Dirac distribution. We know that 5 is the unit operator of 
the convolution. 

Consequently, the response r (M) to signal s(M) (scheme IV. 2) is 

s 

written as r(I) = s(M) ★ h(M) ★ <5(M) , which becomes, for the translated 
signal s(M - Mi): 


. The sum of the signals, SnS^M - J, will 
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s(M-M,) 


r s (M-M- 5 ) 

y 


Ml/ 

s(M) 


r s fM) 


OBJECT PLANE IIo IMAGE PLANE II 


rj-(M) = s(M - Mi)* h(M) ★ 5(M) 
= s(M)* h(M) ★ 6(M - M t ) 
= r g (M - Mi) 

or the translated response. 


Figure 38. Automatic Positioning 
of the Response. In coherent 
optical filtering, for any given filter, 
it is deduced from the properties of 
linearity and invariance by translation 
of the percussional response of the 
optical system that a) a translation 
of the signal in the object plane 
produces the same translation of the 
response in the image plane; and b) 
the sum of n signals gives the sum of 
n responses. 


The physical meaning of this is 

clear. 


The filtering process is written for 
the nontranslated signal s ( M) : 


F T 

a) spectrum analysis, s(M) 1 — ♦ S(£2) ; 

b) filtering by H(Q) , R(O) = S(fi) H(fi) 

c) restitution, R(Q) F ' T * ► s(M)* h(M) 

and for the translated signal s(M) - Mi) : 

F T 

a) analysis, s(M - Mi) ' > S(Q) exp(-2 7ri M t ) ; 

b) filtering, S(£2) H(.S2) exp(-2 7ri O. Mi) ; 


c) restitution, R(£2) exp(-2 7ri Mi) ^r(M-Mi). 

This means that translation Mi introduces, in the Fourier plane, a 
linear phase term containing Mi, which is preserved after filtering and which 
places the response in position. 

Figure 52 gives an example of recognition by a filter designed for 
identifying the character "e": all the "e's" of the text are identified with their 
relative position. 
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3. The Self-Correlation Method 

One method of approaching the problem of pattern recognition 
consists in optimizing the signal/ noise ratio in the exit plane of the system. 
This leads to the theory of the adapted filter [22] . 

However the adapted filter can only operate correctly if the noise is not 
correlated with the signal, or is correlated only to a slight extent. Especially 
for the recognition of letters it is possible to obtain, in the case of similar 
shapes, a response that is more perceptible for the noise than for the signal 
sought. 


It is only in certain limited cases that we obtain an indication of 
presence with a high probability. Nevertheless, it is necessary to discuss 
the method of self-correlation, because the adaptation of this method which 
we have made, notably by constructing hologram filters with low aberrations, 
will be useful in the differentiation method discussed in Paragraph IV. 6. This 
method makes it possible to strongly decrease the probability of false alarm 
when the correlation between the signal and the noise is high. 

The adapted filter . It is shown [16,23] by a reasoning analogous to that 
used in radar technique for the treatment of temporal signals, that the two- 
dimensional filter which optimizes the signal/ noise ratio in the search for a 
signal s(M) is proportional to 

(IV. 4) H(S2) = S* (S2)/ |B(0) | 2 , 

or the quotient of the conjugate complex of the Fourier transform of the signal 
and the spectral density of the noise. 

In Paragraph IV. 5 we shall show that function 1/ |B(S2) I 2 cannot be 
realized physically; hence we are led to the approximation 

(IV. 5) |B(£2) | 2 = constant 

and filter (IV. 4) becomes 

(IV. 6) H(fl)=S*(S2) 

which gives the percussional response 

(IV. 7) H(J2) — — *• g* ( - M) ; 
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hence if, in response r (M) = s(M)^ h(M) due to the signal (scheme IV. 2) , 

s 

we substitute S* (- M) into h(M) , we obtain 

(IV. 8) r (M) = s(M)'* sf* (- M) 

s 

and, in an explicit form, 


which is the self-correlation function of the signal. We know that a self- 
correlation function has a maximum modulus at the origin, or 

r (0) = f \s(m \ I 2 dM . 
s J n v o/ o 

o N ' 

To be sure, the translated signal s(M - Mj) will give a translated signal 

r (M - Mi) having the same maximum but situated at the point M = Mi instead 
s 

of at the origin. 

Likewise, response r^ = s(M) * h(M) due to the noise becomes 
(keeping in mind IV. 7) r fe = b(M) ★ s* (- M) or, in the explicit form, 

r. = f„ b M s* (m - m) dM . 
b -'ll o \ o Jo 
o 

What conclusions can we draw from these expressions? 

a) In the case of a white noise condition (IV. 5) is satisfied and 
the signal/ noise ratio ( IV. 3) becomes [23] 

energy in the signal 

p spectral ( energy) density of noise 


In actual practice it can be said that condition (IV. 5) is assured when 
the signal is only slightly correlated with the noise. Figure 43 gives an example 
of the extraction of words ( signal) which are only slightly correlated to the 
noise (text). 
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b) If the noise is not white, which is often the case in practical 
problems, no conclusion can be drawn at all. It will be seen ( Paragraph IV. 5) 
that the response due to the noise may then be more perceptible than that due 
to signal. This will be prevented in the differentiation method. 


4. The Experimental Arrangement 

The device used both for the differentiation method and for the self- 
correlation method consists a) of a holographic assembly based on [23] , 
which furnishes Fourier filters, and b) a restitutor which identifies the signals 
with the aid of the preceding filters. 


a. The Fourier Holographic Filters 



Let us take a model s 
is to be detected. We have 



of the signal whose presence in object 


s( M o) F ‘ — 


To make the identification we need (expression IV, 6) a filter whose transparency 
is the complex function t = S* (S2) . Except for certain very simple cases, it is 
impossible to realize such a function by currently available means (absorbent 
and dephasing deposits) . By contrast, the holographic technique lends itself 
well to this operation [23] , 

To this end an assembly is constructed (Figure 39) which gives hologram 
filters with very low aberrations. 

A spectrum analyzer, identical to that of the restitutor of Figure 34, 
gives in plane the Fourier spectrum of a signal model s^M q ^ situated in 

plane n^. A part of the laser beam serves for forming a reference wave (shown 

cross-hatched) which is superimposed on the spectrum. The following reason- 
ing is valid for any position of the reference source M , but it will be assumed, 

P 

to simplify the proofs, that point M is situated in the object plane n^. 
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Figure 39. Hologram Filters: Recording Assembly. A spectrum 
analyzer, identical to that of the restitutor of Figure 35, gives in 
plane 5 the Fourier spectrum of the signal situated in object plane 
II Part of the laser beam £ serves for forming a reference beam 
(crosshatched) originating from point Mp situated in plane ll 0 . 

The photographic recording of the phenomenon in plane JF is a 
Fourier hologram. 


1) The signal diffracts on the spectral plane not the Fourier 
transform S(f2) of s(M q ) but (Appendix IV. 1) 

(IV. 11) A q (S2) = B(.0) exp ^i (p\ , 

the phase term exp ( i <? 0 ) is independent of the nature of the signal; it is 

automatically compensated in a restijtutor operating by double diffraction. 
However, in the present case where we use solely the spectrum analyzer part 
of the restitutor, this term persists and must be retained in the formulas. 

/ 

2) The carrier wave originating from point M gives on 

plane (F the following amplitude, standardized to module 1: P 

(IV. 12) Ap(.ft) expj'i <P p j- 

The phase <p (fi) characterizes the spherical reference wave, i. e. , the 
P 

position of point M^. 

3) The photographic plate in fF records the total illumination 
due to (IV. 11) and (IV. 12) . The hologram obtained in this manner presents, 
according to expression (II. 10) , the transparency 


(IV. 13) 


t = 1 - 


Y 


A* A 
P o 


_ JL 


2 


A A 
P 


* 

o 
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which becomes the following expression, keeping in mind ( IV. 11) and (IV. 12) 
with the variable 0 instead of P used in Chapter II, and neglecting the constant 
-y / 2 which is of no importance, 


(IV. 14) 


t * 1 + S(S2) exp 


l i*o 



+ S*(S2) exp 



Let us make 
39). 


<P, 


and w 

P 


explicit by taking S as the origin of the phases ( Figure 


According to Appendix 1, <p Q represents the phase on of a spherical 
wave of radius d having its center at 0, or 

(IV. 15) <p = 7r A d I2 2 . 

Likewise, but by definition, <p represents the phase on of a spherical 

wave originating from point M which -- let us recall — is situated in plane II . 
Whence, by a simple calculation, 

(IV. 16) cp = 7T A d ft 2 - 2 7T S2.M . 

P P 

As a result (IV. 14) becomes 

(IV. 17) t = 1 + S(£2) exp^2 7rifi.M p S* (S2) exp^- 2 rifi.M^. 

Note that, independently of the use as a filter, this hologram makes 
it possible to reconstitute the recorded object or in this case the signal model 
to be identified. 


Figure 4G shows such a restitution for the signal consisting of the word 
"filtrage. " The two images — direct and conjugate — are both virtual because 
of the arrangement of Figure 39. Such a hologram is called a Fourier hologram, 
because the recording is carried out in the Fourier plane. Hence it is necessary 
to use the whole surface of the hologram in order to reconstitute the object 
correctly, contrary to the holograms in diffuse light. In effect, with only one 
portion Afl, only the corresponding frequency band of the object can be recon- 
stituted. 


b. Restitution (Identification of the Signal) 

In the restitutor (Figure 41) the object (= signal + noise) 
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Figure 40, Percussional response of the filter used for the recognition 
of the word ‘filtrage’ in Figure 43. 

( IV. 18) f(M o )= £ (M o ) + b(M o ), 

is introduced in object plane n . 

o 



Figure 41. Pattern Recognition by Self-Correlation. The signal to be 
identified (in the presence of noise) is introduced in object plane ll 0 
of the resti tutor. The Fourier hologram filter is arranged in the spectral 
plane 3\ In the image plane we obtain, as in the case of an ordinary 
hologram, three spatially separated terms: r s (M) = self-correlation 

+ noise; s(M) = signal + noise; r^(M) -- self-convolution + noise. The 
useful term is the self-correlation function of the signal. 


Filter (IV. 17) is introduced into the Fourier plane 3\ 

In that case three responses are obtained in the image plane II, 
corresponding to the three terms of (IV. 17) , and in particular the response 
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sought, i. e. , the sum of the self-correlation function of the signal r (M) 

s 

s(M) ★ s* (- M) and the noise-signal correlation function r, (M) 
b(M)* s*(- M) , b 

To prove this, it is sufficient to follow the filtering scheme (expression 
IV. 2) or, by doing our reasoning on the signal, for example, 

(IV. 19) s(M) —1 , S(fl) ®— S(S2) H(fl). F,T ‘ ?; -»r ( M) . 


H(S2) 

Function H(S2) is the expression (IV. 17) and, since the system is linear, it 
is sufficient to consider the three terms of (IV. 17) separately. 

(1) The Self-Correlation . The third term of (IV.T7) is 

H 3 (I2) = S* (S2) exp(- 2 ir i S2.M p ); 

but 

S(S2) S* (12). F,T< 1 ► s(M) ★ s* (- M) = r (M) 

s 

and 

(IV. 20) S S* exp 2 7ri J2.M —— T — > r (M - M V. 

P s \ p/ 

Thus the contribution of the third term of (IV. 17) is the self-correlation 
function sought but, instead of having for center the origin of plane n, it has its 
center at point 

(IV. 21) M = M , 

P 

i. e. , on the image which would be produced by the reference source on plane 
II if this source were present, , 


On Figure 41 the magnification is not standardized at y = + 1 as in the 

formulas, but where y = - 1, r (M) has its center in M = - M . 

s P 

(2) The Self-Convolution. The second term of (IV. 17) is 


H 2 ( 0) = S*(fi) exp ^2 Tri fi. M ). 
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Let us now set 


(IV. 23) 


rj(M) = s(M) ★ s(M) . 


In plane II we obtain 


(IV. 24) 


S 2 (fi) exp ^2 7 r i Q.M^ F ‘ r t + 



r t ( M + M 


) 


is the self-convolution function of s(M) , with 


its center in 


(IV. 25) M = - M 

P 

which is symmetrical to the center of the self-correlation function. 


The term "having its center" is actually inaccurate because, while a 
self-correlation function exhibits a maximum at the origin (here M j and 
possesses the symmetry property ' 

(IV. 26) |r (m - M Yl = |r (m - m) I, 

s \ P/ s\ p / 

this is not so for the self-convolution function which possesses these two 
properties only for a symmetrical s(M) . 

(3) The Central Term . The first term of (IV. 17) gives 
(IV:27) S(ft) F -— ► s(M) . 


The signal is re-encountered in the center of plane II. In reality, in expression 
(IV. 13) of the transparency of the hologram we have neglected ( Paragraph 
II. 2) certain terms, the result of which being that we obtain a slightly perturbed 
signal. 

The three responses are separated in plane II, provided that we take a 
sufficiently large inclination of the carrier. In actual practice this occurs fo;r 
(Figure 39) 


(IV. 28) OM >2 1, 

P 

where 1 is the half-object field. 


In the same fashion, the noise in the object (IV. 18) will give the three 

terms 
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(IV. 29) 

b(M) 

in the center. 

(IV. 30) 

b(M) ★ s*(m - M ) 

around M , 
P 

(IV. 31) 

b(M) ★ s^M + M p ^ 

around - M . 

P 


c. The Experimental Results 



Figure 42 shows the assembly of the filters and of the restitution 

device. 


Figure 43 gives an example of the recognition of words in a text; the 
responses are the maxima of the corresponding self-correlation functions. 

Figure 44 shows the "noise" due to successive superimpositions of the response 
(recognition of the word "filtrage") . 

It can be seen that for these examples the self-correlation method works 
well, because the words to be recognized are little correlated with the rest 
of the text, but it is quite obvious that if the noise-signal correlation is great, 
there will be ambiguity. This is notably the case with the recognition of isolated 
letters in a text; for example the identification of the letter I by self-correlation 
will also give undesirable responses for the letters which contain the I, i. e. , 

B, D, E, F, H, K, L, M, P, R, T. 

This circumstance is analyzed in Paragraph IV. 5, and the solution which 
avoids the ambiguity is discussed in Paragraph IV. 6, 

d. Notes 


(1) Realization of the Fourier Filters. If during the recording 

of the filter the reference source M is situated in the object plane II ( Figure 

P ^ 

39) , the three responses of Paragraph IV. 4. b will form, on reconstruction, 

in the same plane IL However M may be assigned any given position outside 

P 

plane II o , hence holographic lenses are recorded as a result of which the three 
responses are formed in different planes. 

( 2) The Aberrations of the Filters . If the optical system 
presents aberrations, the signal/ noise ratio maybe considerably reduced. 

The same applies if the Fourier filters are aberrant. 

To study the influence of the filters in this sense we must consider the 
assembly consisting of filter and restitutor. 
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Figure 42. The restitutor and the assembly of Fourier filters. 



Figure 43. Identification of a word in a text by self-correlation: a) 
identification of the word ‘des/ b) identification of the word ‘filtrage. 
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Le filtrage des 
frequences spat i a I es 



Figure 44. Aspect of Noise. Successive superimpositions (from bottom 
to top of the figure) of the response in the identification of the word 
‘filtrage.’ During the measurements, results of the type shown in Figures 
55 and 56 are obtained. 


In the study carried out so far, all the phase terms have been developed 
to within the second order, which amounts to resorting to the paraxial approxi- 
mation. 


Now, it may be shown that, without carrying out any approximations and 

for any given position of the reference source M , we obtain the useful response 

P 

(IV. 20) with highly reduced aberrations, within the limits of validity of the 
Fourier formalism in optics. 

To this end it is necessary that this formalism be effectively realized: 
this is the case of the device with convergent light which has been used 
(Appendix 2) . Moreover it is necessary, for the correction of the holographic 
aberrations, that the spectrum analyzers of the restitutors and of the recording 
device of the filters be identical. In particular , the coma and the spherical 
aberration are then strictly zero. 

(3) The Formalism in Pattern Recognition . In this case, it 
is sufficient, on the basis of the preceding results, to consider only the useful 
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response (IV. 20) and the corresponding noise (IV. 30) , and to neglect the 
translation (due to the inclined reference wave) of the response in the image 
plane. The filtering scheme is then written as 


s(M) + b(M) — — — — «-S(£2) + B(£2) — ® F ' T ~ 1 
signal noise spectrum j 


r (M) 
s 

+ r, (M) 
b 


= s(M) * s* (- M) 
= b(M) * s* (- M). 


(IV. 32) 


S* (£2) 
filter 


5. Introduction to Identification by Differentiation 

To show how one may identify a signal with a low probability of 
false alarm, it is necessary first of all to examine the reasons which bring 
about failures in the classical self-correlation method. 

a. The Adapted Filter 


Let us take an object f(M) = s(M) + b(M) consisting of a signal 
s(M) to be identified — e. g. , a printed letter — in the presence of what, by 
extension, will be called noise b(M) , i. e. , all the elements of the object other 
than s(M) (the rest of the text) . 

When we identify the signal by means of a filter of transparency S* (£2) , 
as in (IV. 32) , we obtain in the image plane II a self-correlation function r 

s 

whose maximum ought to indicate a) the presence of the signal, b) the position 
of the signal ( Cartesian coordinates of the maximum) . 

However, we obtain also the functions of the signal/ noise intercorrelation 
r^, and the latter may be much more perceptible than r g . Thus a number of 

false responses are obtained in this manner. 

hi effect, in Paragraph IV. 3 we have started out from the adapted filter 

(IV. 33) H(£2) = S* (£2)/ |B(£2) I 2 . 

This filter gives good results in radar technique because the noise, 
which originates essentially from the reception chain, is known and it can often 
be considered as white in the reception band. 
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By contrast, in pattern recognition and, in particular, in the case of a 
printed text, the conditions are quite different. While the numerator of (IV. 33) 
may be realized by holography (Paragraph IV. 4. a) , the denominator 1/ |B(fi) | 2 
cannot generally be realized. In effect, 

1) The noise b(M) is unknown, because we consider as noise 
all the unwanted letters of the text, and this is unknown; 

2) The noise actually consists of undesirable signals of the 
same type as the signal to be identified; they have a limited support and their 
spectrum B(Q) presents numerous zeros. Assuming the noise to be known, 
function 1/ |B I 2 can be realized only in a very approximate manner, by replacing 
an infinite transparency by a unit transparency ; 

3) even if we assume the noise to be known, it would be 
necessary to construct an antinoise filter for each page of the text. 

Under these conditions we are led to the approximation 
(IV. 34) |B (fi) I 2 = constant 


Filter (IV. 33) becomes H(J2) = S* ( fi ) » and the response becomes, 
according to (IV. 32) : 

(IV. 35) r(M) = r (M) + r (M) , 

S D 

whose first term 

(IV. 36) r (M) = s(M) * s* (- M) , 

s 

due to the signal is the self-correlation function of the signal, and the second 
term 

(IV. 37) r b (M) =b(M) * s* (- M) 

is the signal/ noise correlation function. 

b. The Causes of Failure. 

When the noise is highly correlated with the signal, the identification 
may become difficult and at times even impossible. 

In effect, expression (IV. 37) may be more readily perceptible than the 
self- correlation function (IV. 36) . 
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This is particularly striking in the extreme case where the noise con- 
tains the signal (example: the letter E contains the letter I) , For simplifica- 
tion we shall assume that the signal and the noise have an amplitude of 1 in the 
regions s(M) and b(M) of respective areas 2 g and 2^ , with 2 g < 2^. 

Furthermore, let us take c(M) as a region common to s(M) and b(M) . 

To sum up, then, 

(IV. 38) s(M) C b(m) , c(M) = s(M) H b(M) . 

From this it can be inferred that r g and r^ have the same maximum, which is 
precisely 

(IV: 39) / ls(M) I 2 dM. 

o 

However, while r. (M) is maximum in a region having an area 
b 

(IV. 40) 2 = 2- 2 , 

b s 

the self-correlation function r (M) is maximum only at point M = 0. This 

s 

means that while we look for s(M) we identify b(M) , 
c. Example 


Let us illustrate this by an example. Take a text consisting of the 
two letters I and E ( Figure 45a) . 

To simplify this demonstration we shall carry out our reasoning in one 
dimension — x — only. 

Let the object be 


(IV. 41) f(x) = s(x) +b(x) 

the signal to be identified is 
(IV. 42) I-s(x) 

and the noise 

(IV. 43) E = b(x) . 
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The distribution of the amplitudes in the object may be represented, 
along x, by two gates of different widths (Figure 45b) . 


SIGNAL 


I ' 
SIGNAL 

s(x) 

I 



RESPONSE X 


s(x)*s*( -x) 



NOISE 


NOISE x 

b(x) 

I 


(b) 



Figure 45. Disadvantages of the Method of Pattern Recognition by Self- 
Correlation. The object (a) gives the distribution of amplitudes (b) along 
x. If the object of line (b) is correlated with s(x), the response (c) consists 
of two terms: 1) a triangle which is a self-correlation function of the signal, 
and 2) a trapeze, which is a function of the signal/noise correlation. 

The second term, which is much more perceptible than the first, has the 
result that it is the letter E that is identified while the letter I is sought. 
Compare this result with that of Figure 46. 


The response (Figure 45c) consists of a triangle, a self-correlation 
function of the signal, and a trapeze, which is the signal/ noise correlation 
function. 

We note that the trapeze is much more readily perceptible than the 
triangle which, in principle, ought to permit the identification of the signal. 
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6. The Differentiation Method 

a. Principle 

Expression (IV. 40) shows that the ambiguity in the recognition 
is due to the finite surface of the elements considered. Now, a shape is 
generally defined by its contour and, more particularly, by its lines of dis- 
continuity. The zones of constant or slowly varying amplitude play a minor 
role if it is desired to only recognize the presence and position of an object; 
thus they may be considered as a "noise" which one may get rid of, notably 
by differentiation. 

Hence a suitable differential operator is applied to the object and the 
differentiated object is correlated with the signal to which the same operator 
has been applied. The differentiation is obtained simply by placing in the 
spectral plane fF an absorption filter having an appropriate law of transparency. 

This will be made more explicit below in the one-dimensional and in 
the two-dimensional case. 

For the following discussion let us remember that the experimental 
device used is the same as in Paragraph IV. 4, and that the filtering scheme 
adhered to is still (IV. 2) . 

b. The One-Dimensional Case 

The principle of the method will seem more obvious if we went 
back to the example of Paragraph IV. 5. c relating to the case of simple self- 
correlation. 

The object ( Figure 46a) is differentiated ( Figure 46b) . 

(IV. 41) f'(x) = s'(x) +b'(x), 

then it is correlated with the differentiated function s Pli (x) , which gives the 
response 

(IV. 42) r(x) =r g (x) + r fe (x) 

consisting of the response due to the signal sought 

(IV. 43) r (x) = s'(x) * s'* (- x) 

s 

and of the response due to the noise 

(rv. 44) r fo (x) = b'(x) * s'* (- x) . 
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SIGNAL NO'SE 


X 

DIFFERENTIATED SIGNAL 
s’(x) 


x 


RESPONSE 

s’(x)*s’*(-x) 


Figure 46. Pattern Recognition by Differentiation. The object (a) shown 
in the figure is the same as that of Figure 45. The differentiated object 
is shown in (b). If this object is correlated with a signal model which 
itself is differentiated, we obtain (c) where the useful response exhibits, 
with respect to the noise, a signal/noise ratio of 2 in amplitude and 4 in 
intensity. 

The real signals are generally two-dimensional, in which case the signal/ 
noise ratio is 16 for a similar example (Figure 47). 

We then see on Figure 46c that the amplitude of the maximum of r is 

s 

double its secondary maxima, and likewise the double of all the maxima due to 
the noise. 

Hence the signal/ noise ratio is 4 in intensity. 

The objects generally considered are two-dimensional objects. For a 
similar example in x, y, the signal/ noise ratio is then 16, instead of 1 as in 
the self-correlation example of Paragraph IV. 5. c. 

c. Form of the Fourier Filter 

To determine the nature of the filter which leads to the results of 
the preceding paragraph, we must consider the Fourier plane. Let u be the 
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reduced variable in this plane. We then have, starting from the desired 
response, 

(IV. 45) s’(x) s'* (- x) — F 1 — 4 n 2 u 2 S(u) S* (u) , 

(IV. 46) b’(x) * s'* (- x)- F ‘ 4 Ti 2 u 2 B(u) S* (u) ; 

hence the filtering scheme (IV. 2) becomes ( except for the coefficient - 4 71 2 ) 

(IV. 47) f (x) — L f(u)- ® -* ->r(x) = f'(x) * s'* (- x) 

u 2 S* (u) 

which shows that the filter consists of the product of u 2 and S* (u) ; i.e, , 
brought about materially through juxtaposition 

a) of a filter of transparency S* (u) , which may be realized by 
holography ( Paragraph IV. 4. a) , and 

b) a filter whose law of transparency in amplitude is 

(IV. 48) t = u 2 , 

or a parabolic transparency having a minimum 0 in the center. 

This is possible, because the spectral plane O' is limited by the entry 

diaphragm of lens L. If u is the radius of this diaphragm, the filter has the 

m 

transparency in amplitude 

(IV. 49) («/u m )» in the diaphragm. 

0 outside the diaphragm. 

d. The Two-Dimensional Case 

To transpose the preceding reasoning to the two-dimensional case, 
it is necessary to differentiate the object f(M) and the pereussional response 
s* (- M) of filter S* (£2) with respect to x and y . Then any differentiation will 
give a result. However, it is clear that the optimum operation is that which 
makes the "lines of the greatest slope" appear in the luminous distributions 
f(M) and s* (- M) . This is obtained by the gradient operator 

(IV. 50) f (M) =>grad*f(M) , s* (- M) =>gracT s* (- M) . 
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This gradient apparently has no physical meaning if considered in 
isolation. However, this is not so if we consider the whole recognition process 
in its entirety. 

It is then shown (Appendix 3) that, if the diaphragm of lens L of the 
restitutor has a radius 0 , the two-dimensional filter has the form (also 

Paragraph IV. 6. c) m 

jjy 51) H(S2) = (fi/fi m ) 2 S* (fi) in the diaphragm 

H(S2) = 0 outside the diaphragm; 

the response due to the signal then becomes 

(IV. 52) r = / grad s(M ) . grad's* (M - M) dM 

o J J 

and the noise gives 

(IV. 53) r, = L grad b(M ) . grad s* (M - M) dM . 

b J a o o o 

o 

These expressions are independent of the coordinate axes, and the 
detection of the lines of discontinuity depends neither on the form nor on the 
position of the objects, because the expression under the integration sign is a 
scalar product of gradients. 

e. Differentiations of a Higher Order 

In the same manner, if necessary by using a filter with (£2/S2^) 4 
we obtain for the signal (A representing the Laplacian operator) 

(IV. 54) r = As(M) As*(- M) 

s 

and for the noise 

(IV. 55) r, = Ab(M) Jf As* (- M) 

b 

The order of differentiations is limited only by the power of the lasers employed, 
because the differentiating filters have a maximum absorption at the center. 

NOTE: An approximate differentiation, which consists in carrying out 
a high-pass filtering, gives satisfactory results in certain cases. 
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f. Experimental Results 


Using the experimental device of Paragraph IV. 4 and the filters of 
form (IV. 51) , the following results have been obtained: 

Figure 47 is the two-dimensional illustration of the gate pulses of 
Paragraph IV. 6. b: a) the signal to be identified is indicated by an arrow; 
b) the response shows the presence and position of the signal without ambiguity. 
The signal/ noise ratio is 26 dB (on the recording made with the aid of a film 
of y ^ 2) . 


Figure 47. Method of Pattern 
Recognition by Differentiation 
The signal to be identified in 
object (a) is indicated by an 
arrow. In the response (b) the 
signal/noise ratio is 26 dB. 


Figure 48 shows two cases of 
fingerprint identification: (a) and (b) 

showing the signal, and (c) and (d) the 
c orresponding responses . 

Figure 49 relates to the identifi- 
cation of fingerprint (a) in the presence 
of print (b) having a very similar form; 
(c) is the response. 


Figure 50 shows a fingerprint 
and its spectrum; we can see that it is 
difficult to realize such a Fourier 
spectrum by a method other than 
holography. 

Figure 51 shows the identifica- 
tion of a letter in a text: in text (a) the 
signal to be identified is the letter "e"; 
the response is shown in (b) , where all 
e's are simultaneously identified accord- 
ing to the properties of Paragraph IV. 2. b 
(automatic position-locating of the 
response) . 


Figure 52 shows an example 
analogous to that of Figure 51, but with 
different characters; in object (a) the 
signal is letter "e, " while (b) is the 
response. We shall return to Figures 51 
and 52 later on. 

Figure 53 likewise shows the 
extraction of letter "s." 
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Finally, Figure 54 shows the extraction of letter O in the presence of 
two similar letters, C and Q; the response in (b) indicates a signal/ noise ratio 
of 4 dB. 


g. The Signal/Noise Ratio 

This is greatly improved in the differentiation method, as we 
have seen. However this ratio always remains a function of the experimental 
conditions, and especially of the aberrations. 

Figure 55 shows a microdensitogram of a line of the response (b) of 
Figure 51. 

The text is plotted at the bottom of the figure, the letters being 
arranged opposite the corresponding peaks of the response; it is noted that the 
signal/ noise ratio expressed in difference of density between the minimum "e" 
signal and the maximum noise is 0. 12, or 1.2 dB. 

The experimental device was equipped with good commercial lenses. 

Figure 56 represents an analogous microdensitogram for the response 
of Figure 52; the signal/ noise ratio has risen to 8 dB. The principal reason 
of the improvement is the use of a new optical system (calculated by the CERCO 
Company) which has been well corrected for aberrations; the latter are shown 
on Figure 57. 
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Figure 48. Differentiation Method. Identification of a fingerprint. 
Recognition of fingerprints (a) and (b) in the presence of fingerprints 
(a’) and (b* ). Responses (c) and (d) indicate the presence and position 
of fingerprints (a) and (b). 



Figure 49. Differentiation Method. 
Identification of a fingerprint (a) in 
the presence of a very similar finger- 
print (b). The response is shown in 
(c). 


7. Conclusions on Pattern 
Recognition 

Pattern recognition 
by differentiation consists merely 
in the juxtaposition, in the spec- 
tral plane of a restitutor oper- 
ating by double diffraction, of 

a) A filter whose trans- 
parency is proportional to the 
conjugate complex of the Fourier 
transform of the signal to be 
identified; this filter may be 
.realized by holography; 

b) An absorbent filter 
whose law of transmission 
entails the application of a 
differential operator to the object 
( signal + noise) and to the per- 
cussional response of the pre- 
ceding filter. 
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Figure 50 o Example of an Object and of its Fourier Spectrum: a) object; 
b) spectrum. 


le filtrage des frequences spatiales 
le filtrage des frequences spatiales 
le filtrage des frequences spatiales 
le filtrage des frequences spatiales 
le filtrage des frequences spatiales 
le filtrage des frequences spatiales 



Figure 51. Method of Recognition by Differentiation. In object (a) the signal 
to be identified is the letter s e/ The response (b) shows the simultaneous 
identification of all the Vs* of the text (linearity of the optical system): the 
responses are formed automatically at the site of the corresponding letter. 
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le filt rage des frequences spatiales 
le filtrage des frequences spatiales 
le filtrage des frequences spatiales 
le filtrage des frequences spatiales 
le filtrage des frequences spatiales 
le filtrage des frequences spatiales 



Figure 52. Method of Recognition by Differentiation. The signal to be identified 
in object (a) is the letter ‘e.’ Response shown in (b). 
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le fsltrage des frequences s pa h a les 
le fill rage des frequences s pah ales 
le filt rage des frequences spatiales 
le fill rage des frequences spatiales 
le fi If rage des frequences spatiales 



Figure 53. Method of Recognition by Differentiation. The signal to be identified 
is the letter ‘s’ in text (a). Response is shown in (b). 


If necessary, these two filters may be realized on a single support. 

Efforts in this area have been directed essentially at the identification 
of characters in a printed text. The results obtained in this difficult case 
(the letters have very similar shapes and a character represents a small part 
of the total energy of the text) make it possible to envision the extension of the 
method to other types of signals, i. e. , electrical or acoustic signals. 
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COQ 



Figure 54. Differentiation Method. Identification of the 
the letter 0 among two other, very similar letters: a) 
text; b) response. 



le filtrage des frequences spatiales 


Figure 55. The Signal/Noise Ratio.- On this microdensitogram of a line 
of the response of Figure 51, the text is plotted at the bottom of the figure, 
and the letters to be identified are underlined. The optical densities are 
plotted along the ordinate axis; The rather mediocre signal/noise ratio is 
due to the aberrations of the optical system and of the hologram filters. 
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lefiltrage des frequences spatiales 


Figure 56. The Signal/Noise Ratio. This example, analogous to that of 
Figure 55, relates to the extraction of letter ‘e’ in Figure 52. The improve 
ment of the signal/noise ratio is due to the correction of the aberrations of 
the optical system and of the hologram filters. 
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APPENDIX 1 

Analysis of the Filtering of Spatial Frequencies 


It is attempted to carry out optically the operation whose principle is 
schematically shown below: 



filter 


by means of the arrangement represented in Figure 58. 



It is shown [19] that to ensure linear filtering, the restitution lens L 
must have its entry diaphragm in the spectral plane This is schematically 
represented by a simple lens placed in this plane. 

1. Spectrum Analysis . Under the effect of the convergent coherent 
illumination originating from lens L q , the signal f * n Pl ane n 0 diffracts, 

on plane & having a general point P(x ls y*) , the amplitude (expression 1. 16) 

( A . 2) g( P) = exp ( i k P 1 2 /2 d) F ( P) , 

where d is the distance between planes and ff, and F(P) the Fourier 
transform of f . In reduced coordinates, with 

(A. 3) .0= P/Xd, 
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expression (A. 2) becomes, by writing g(0) and F(0) for g(A d 0) and F(A d 0) , 
(A. 4) g(0) = exp(i 7 r A d 0 2 ) F{0) ; 

this spectrum is the product of the Fourier spectrum 

(A. 5) F(0) = f n f exp^- 2 7r i 0. dM Q 

and the quadratic phase factor 

(A. 6) exp(i ir A d 0 2 ) . 

2. The Filtering. After passing through the filter of transparency H(0) 
placed in plane , wave (A. 4) becomes 

(A. 7) exp(i 7r\d fi 2 ) F(fi) H(S2); 

hence term (A. 6) has no effect on the filtering. However it must be compensated 
in order that the restitution operation give a good inverse Fourier transform of 
the product F(fi) H(0) , 


3. The Restitution . It is known that a lens of focal length 1 introduces 
in its plane (plane 5) the phase difference 

(A. 8) exp(- i k P 2 /2 1), 

or, in reduced coordinates 

(A. 9) exp(- i 7r A d 2 0 2 /l) . 


After passing through lens L, the wave (A. 7) consequently becomes 


(A. 10) 


gi(0) = exp 


i 7r A d 2 0 2 



F(0) H(0). 


This wave gives, by diffraction on a plane at a distance d' from 3 r , an amplitude 
distribution r(M) , which is obtained by application of formula (1. 7) relating to 
the Fresnel diffraction. This is written, in reduced coordinates, as 



whence, by expanding the exponent and taking (A. 10) into account, 


(A. 12) r(M) = exp (i j f e icf> F(fi) H(fi) exp 


-2 7r i — O.M Idft, 


with 


(A. 13) 


?=nd2fi2 (T + i-Tj • 


Now, since 11^ and Ilare conjugate planes, we have 
( A. 14) “d + "d 7 " T = ° ; 


hence 


(A. 15) e 1 ^ = 1. 

Moreover, (A. 12) is the final result of the filtering; at this stage the recording 
is made by a quadratic detector and the phase term 

(A. 16) exp(i 7 r M 2 /\ d’) 

disappears. It may therefore be neglected. Hence (A. 12) becomes 

(A. 17) r(M ) = / F(O) H(S2) exp (- 2 ni-^ 

which represents the inverse Fourier transform of the product 
(A. 18) R(O) - F(O) H(S2) 

which, except for a scale term, is magnification \ = - d'/d of lens L. 
Standardizing y at 1, we have, as a result, 

(A. 19) R(fl) F — 1 — »r(M) = f(M) 4 h(M) . 

Therefore the filtering follows the schematic outline of (A. 1) , 
with the phase term (A. 6 ) having been compensated and term (A. 16) dis- 
appearing during the recording. 


86 



Nevertheless, in Fourier transform holography it is necessary to keep 
(A. 6) in mind, and if response (A. 12) must be taken up by another coherent 
system, it is also necessary to take (A. 16) into account. 
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APPENDIX 2 

Definitions Relating to the Photographic Emulsions 


We define the following parameters for the photographic emulsions: 

1. Transparency in Intensity . This is the ratio, T, of the emergent 
fluxes and the incident fluxes $5 at each point of the emulsion on which the 
image is prepared: 

(A. 20) T=V$o. 0STS1, 

Thus the transparency T varies with each point of coordinates x, y on 
the emulsion, and it defines a function T(x, y) of these coordinates. 

2. Optical Density . 

(A. 21) D = log(l/T). 

3. Lumination. This is the energy W received by the plate, or the 
product of illumination E and the exposure time t 

(A. 22 ) W=Et. 


4. Characteristic (or Blackening) Curve. This curve gives the varia 
tions of density as a function of lumination ( Figure 59) . This curve often 
includes a straight-line portion whose slope y is also called contrast factor. 
The straight-line portion has the form 


(A. 23) D = y log( W/W Q ) , 

where W 0 is a constant. 



According to (A. 21 ) we then 
have T = (W/W 0 )" Y , or 

(A. 24) T = CE i '; 

E is the illumination received by the 
plate and C a constant which depends 
notably on the exposure time, on the 
sensitivity of the plate and on factor y 


Figure 59. 
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5. Transparency in Amplitude . In coherent optics it is the (complex) 
amplitudes of the waves which intervene; then the transparency in amplitude is 
defined as the ratio (t) of the emergent amplitude (A t ) and incident amplitude 
(A 0 ) at each point of the emulsion on which the image is produced. 

(A. 25) t= Ai/A 0 , 0 £ It I =£ 1. 

t(x,y) is a real or complex function of coordinates x, y on the emulsion; its 
modulus is less than 1. 

We also have the relations 

(A. 26) T = 1 1* 

i 

and, for a real t, t = T 2 , or 

(A. 27) t = C' e”"^ 2 , 

where C' is a constant. 
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APPENDIX 3 

The Two- Dimensional Differentiating Filters 


Let us take an identification by classical self-correlation 
(A. 28) r(M) * f(M) s* (- M) . 

The detection of the lines of discontinuity ( Paragraph IV. 6. 4) leads to the 
determination of the linear differential operators L which, when applied to the 
object 

(A. 29) f(M) m s(M) +b(M) 

and to the percussional response of the Fourier filter 


(A. 30) 


s* (- M) 


give a result 


(A. 31) 


a) independent of the coordinate axes ; 

b) independent of the shape of the signals treated, i. e. , of their 
i position and orientation. 


Now let us express the Dirac distribution by 5(M) = <5(x, y) (see, for 
example [17] for the mathematical formalism) , and let the linear differential 
operator be 


(A. 32) L(M) = L(x,y) = £ (a 


8 m <5 9 n 5 

+ b — — + c 
m _ m n „ n p 


9 P 5 


9x q 9y p - q 


9x 9y 

which, when applied to (A. 29) and (A. 30), transforms (A. 28) into 


(A. 33) 


Vi ( M) = [L*f<M)J * [L*s* (- M) ] 
= L ★ L Arr(M) . 


The Fourier transform in S2(u,v) of r(M) is 
(A. 34) R(S2) = F(fi) S* ( S2) , 
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The Fourier transform of L is 


(A. 35) 


L = 2 


a (27riu) m + b (27riv) n +e (27ri) P u <3 v P P 
m n p 


and that of (A. 33) 


(A. 36) R 1 (Q) = L 2 R (0) , 

To satisfy (A. 31) it is necessary that L 2 have a radial symmetry 

L 2 = a( 2 7r i) 2n (u 2 + v 2 ) n = a(2 7r 1 0) 2n , 

where n is a whole number and a a constant which will be subsequently 
neglected. 

Then the response (A. 33) becomes 

(A. 37) r t (M) = r(M)'*F.T. -1 {- 4 7i*(u s + v 2 ) 

Then it is necessary to distinguish between several cases, depending on the 
value of n. 


1. Example of n = 1. We have 


(A. 38) 


. «/ 2 2\ E. X" 1 • 3 2 6 8 2 5 

- 4 ir (u z + v 2 ) 1 + 7TT • 


3x 3y 

As a result the response (A. 37) becomes (keeping A. 28 in mind) 

Vd 3 2 6 N 


(A. 39) 


Ti(M) = 


dx i 3y 2 


★ f(M) ★ s* (- M). 


This convolution product is distributive relative to the addition and 
to the differentiation, whence 


(A. 40) 


r t (M) = 


8f(M) 

3x 


3s* (- M) 
3x 


3f . 3s* (- M) 

3y 3y 


Now df/dx, df/dy and 3s* (- M)/3x, 3s* (- M)/3y are the components of 

grad f(M) and grad s* (- M) ; hence (A. 40) is written — by making the con- 
volution products explicit — as the sum of a scalar product 
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(A. 41) 


rj(M) = f grad f(M ) . grad s* (M - M) dM . 

11 o o o 

o 

2. Example of n = 2 . We have 

16 ^(u 2 + v 2 ) 2 — T ’ — AS if AS. 

Then the response (A. 37) is written as 

rj (M) = AS * f(M) * AS * s* ( - M) 

or also as 

r^M) = Af(M) if As* (- M) . 
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